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Abstract 

A review of the theoretical and experimental status of hybrid hadrons is presented. The states 
7ri(1400), 7ri(1600), and 7ri(2015) are thoroughly reviewed, along with experimental results from 
GAMS, VES, Obelix, COMPASS, KEK, CLEO, Crystal Barrel, CLAS, and BNL. Theoretical 
lattice results on the gluelump spectrum, adiabatic potentials, heavy and light hybrids, and 
transition matrix elements are discussed. These are compared with bag, string, flux tube, and 
constituent gluon models. Strong and electromagnetic decay models are described and compared 
to lattice gauge theory results. We conclude that while good evidence for the existence of a 
light isovector exotic meson exists, its confirmation as a hybrid meson awaits discovery of its 
iso-partners. We also conclude that lattice gauge theory rules out a number of hybrid models 
and provides a reference to judge the success of others. 


1 



Contents 


1 Introduction |4] 

2 Lattice Gauge Theory [5] 

2.1 Adiabatic Gluonic Surfaces. [5] 

2.2 Gluelumps. [7| 

2.3 Heavy Hybrids. [8] 

2.4 Light Hybrids. [12] 

2.5 Hybrid Baryons. [TH 

2.6 Transitions . M 

2.6.1 Vector Hybrid Mixing. [15] 

2.6.2 Heavy Hybrid Hadronic Transitions . [16] 

2.6.3 Light Hybrid Hadronic Transitions. [IT] 

2.6.4 Hybrid Radiative Transitions. [IT] 

3 Hybrid Models [18] 

3.1 Bag Models. [18] 

3.2 Relativistic String Models. [20] 

3.3 Flux Tube Model. [21] 

3.3.1 Model Development. [21] 

3.3.2 Further Developments. [23] 

3.4 Hybrid Baryons in the Flux Tube Model. [25] 

3.5 Constituent Glue Models . [25] 

3.6 Sum Rule and Schwinger-Dyson Calculations. [27] 

3.7 Model Comparison. [28] 

4 Hybrid Decay Models [29] 

4.1 Flux Tube Decay Model. [29] 

4.2 The PSS Hybrid Decay Model. M 

4.3 Hybrid Photocoupling and Production. [3l] 

5 Experimental Situation [32] 

5.1 Exotic Quantum Number states. [32] 

5.2 The TTi(1400). [33] 

5.2.1 GAMS results on the 7ri( 1400). [33] 

5.2.2 KEK results on the 7ri( 1400). [33] 

5.2.3 VES results on the 7ri(1400). [33] 

5.2.4 E852 results on the 7ri( 1400) . [31] 

5.2.5 Crystal Barrel results on the 7ri( 1400). [33 

5.2.6 OBELIX results on the 7ri( 1400). [33 

5.2.7 COMPASS results on the TTi (1400). [33 

5.2.8 Interpretation of the 7ri( 1400). [33 

5.3 The TTi (1600). M 

5.3.1 VES Results on the 7ri(1600). [33 

5.3.2 E852 Results on the 7ri(1600). [33 


2 





































5.3.3 Crystal Barrel results on the 7ri(1600) 

5.3.4 CLEO-c results on the 7ri( 1600) . . . , 

5.3.5 CLAS results on the 7ri(1600). 

5.3.6 COMPASS results on the TTi (1600) . , 

5.3.7 Summary of 7ri(1600) results. 

5.4 The TTi (2015). 

5.4.1 E852 Results on the 7ri(2015). 

5.4.2 Interpretation of the TTi (2015). 

6 Outlook [53] 


3 












1 Introduction 


For forty years Quantum Chromodynamics (QCD) has served as one of the pillars of the Standard 
Model. Early doubts about the existence of quarks were dispelled by the discovery of the J/'0, 
while—following a suggestion due to Ellis, Gaillard, and Ross - gluons were confirmed by the 
discovery of three jet events at DESY in 1979 [^. 

Although gluons are now firmly established as the carriers of the strong force, their nonpertur- 
bative behavior remains enigmatic. This unfortunate circumstance is chiefly due to two features of 
QCD: the theory is notoriously difficult to work with in the nonperturbative regime, and experi¬ 
mental manifestations of glue tend to be hidden in the spectrum and dynamics of hadrons. Hadrons 
that carry valence quark and gluonic degrees of freedom are one such experimental manifestation 
that has been postulated since the early days of QCD [^. These states are called hybrids and a 
review of their experimental and theoretical status is the purpose of this article. 

That discovering explicit nonperturbative glue is difficult can be gleaned from the history of the 
development the quark model and QCD. All the well-established mesons have equal to 0 
0^^, 1^^, etc, which can be created with fermion-antifermion pairs in an given orbital momentum 
state. This strongly suggests that quarks are spin-1/2, while the spectrum ordering suggests that 
energy increases with orbital angular momentum. Furthermore, the absence of mesons with isospin 
or strangeness greater than unity supports the idea that mesons are fermion-antifermion states. 
Thus the simple quark model of mesons (and baryons) was partly motivated by the absence of 
“exotic” hadrons such as multiquark or gluonic states. It is thus perhaps no surprise that QCD 
exotics have been difficult to observe. 

The simplest explanation for this observation is that gluonic degrees of freedom are somewhat 
“stiff’ and therefore difficult to excite. Of course, with the increasing energy, luminosity, and 
capabilities of modern accelerators and detectors one might hope that this regime can be explored. 
The interesting question then is: can an analogous spectroscopy of glue be discerned in the data? 

A basic question that the data and its interpretation must address is what is meant by the notion 
of valence glue. As an example of the importance of choosing correct degrees of freedom, we mention 
the simple problem of determining the number of components of a constituent gluon. For example, 
one might suppose that a massive constituent gluon (of mass nig) should be transverse so as to 
maintain consistency with Yang’s theorem. However, this is inconsistent with the requirements of 
Lorentz invariance. Thus, for example, J — 1 glueballs are expected to exist and lattice calculations 
indicate that they are quite heavy (roughly 3 GeV) [^. Such a state may not be constructed from 
two transverse constituent gluons (due to Yang’s theorem) and therefore may be expected to have 
a mass of roughly 3m^. However massive vector gluons have no such constraint and one therefore 
expects them to have a mass of approximately 2m^. Clearly the assumed active degrees of freedom 
play a large role in determining the gross features of the spectrum. 

There are two broad ideas concerning soft glue: it is some sort of string or flux tube or it is 
an effective constituent confined by a bag or potential. In different language, nonperturbative glue 
can be thought of as collective, nonlocal degrees of freedom, or as a local quasiparticle degree of 
freedom. Lattice gauge theory and experiment are the two main methods which permit resolving 
this longstanding issue. Of course the primary task is to experimentally verify that such excitations 
of QCD exist. Then a thorough investigation of the hybrid spectrum and its production and decay 
mechanisms must ensue; this must be compared to lattice gauge and model computations, and a 
qualitative picture must be built before we can claim to fully understand the enigmatic gluonic 
sector of the Standard Model. 
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We will start with a review of properties of hybrids and related gluonic systems as determined 
by lattice gauge computations. This will serve to set terminology and establish a base to which 
models can be compared. This is followed by brief summaries of the ideas and main results in 
string-based models, bag models, and constituent glue models. We then summarize the current 
status of the experimental search for hybrid mesons. Finally, future experimental prospects and an 
outlook on what is required and expected from theory is presented. 

Before starting, it will be useful to note that valence gluonic degrees of freedom increase the 
quantum numbers that are available to fermion-antifermion systems. The parity and charge conju¬ 
gation that are available to qq systems are specified by P = C = where L and 

S are the qq total angular momentum and spin. This means that — 0 , odd and even^ 

are not available to simple quark-antiquark systems. Hadrons with these quantum numbers are 
called “(quantum number) exotic”. A variety of notations have developed over the years for hy¬ 
brids. Terminology appears to have settled on using standard PDG notation to describe the flavor 
of a hybrid, a subscript to denote the total spin, with parity and charge conjugation following that 
of the named hadron. Thus, for example, a = 1 ^ isovector hybrid would be denoted tti. 

2 Lattice Gauge Theory 

The use of a lattice regulator for Euclidean field theories permits the numerical evaluation of 
path integrals. Thus all correlation functions can, in principle, be computed under controlled 
approximations. Historically, this brilliant promise has been compromised by various technical 
issues including large quark masses, the difficulty in computing “hairpin” fermion lines, difficulty 
in computing the excited spectrum, and most prominently, difficulty in computing the determinant 
of the Dirac operator. The latter issue was handled by ignoring it (this is called the “quenched 
approximation”; it is, in fact, not an approximation, but a truncation that renders a quantum 
field theory inconsistent). The root problem of all of these issues is that they introduce excessive 
noise into observables - however, to some extent all of them have been overcome in the past 
decade. Thus modern lattice computations are made on large lattices, are not quenched (they 
include the effects of virtual quarks), incorporate large bases of interpolating fields (which permits 
the extraction of excited states), and employ sophisticated stochastic methods to compute hairpin 
fermion propagators. The remaining issues are chiefly quark masses that tend to be too large 
and the necessity of incorporating continuum interpolating operators. Both of these are rapidly 
becoming things of the past. 


2.1 Adiabatic Gluonic Surfaces 

Lattice investigations of soft gluonic matter date to an investigation of the adiabatic gluonic spec¬ 
trum by Griffiths, Michael, and Rakow [^. The idea is that gluons are a “fast” degree of freedom 
with respect to heavy quarks. In the static limit the quark and antiquark serve as a color source and 
sink at a distance R and the gluonic field arranges itself into configurations described by various 
quantum numbers. These quantum numbers match those of diatomic molecules: the projection of 
the gluonic angular momentum onto the qq axis (denoted A), the product of gluonic parity and 
charge conjugation 77 = {PC)g^ and T-parity, which represents the reflection of the system through 
a plane containing the qq axis. The notation 


5 


will be used in the following, where A = 0 , 1 , 2 ,... is denoted S, E, A, ...; 77 is labelled u (negative) 
or g (positive); and y = zb. For A ^ 0 the two F-parity eigenstates are degenerate and the Y 
quantum number is not given. 

Griffiths et al sought to compute the adiabatic configuration energies, E{A^;R) as a function 
of the distance between the color source and sink. This computation was made in SU(2) quenched 
gauge theory on an 8 ^ lattice and rough results for the lowest two surfaces were obtained. 

Subsequent work was carried out by the UKQCD collaboration and Juge, Kuti, and 

Morningstar [^. Juge et al. performed a large scale computation with many basis states and were 
able to obtain many excited states out to large distance. The results are presented in Fig. here 
ro is approximately 0.5 fm and the potentials have been renormalized by subtracting the value of 
F;(y+) at ro. 

An interesting feature of the spectrum is that the adiabatic surfaces tend to finite values at 
small interquark separation. This may be a surprise because the expected perturbative behavior 
at small R for a pair of quarks in a color octet is 


Alternatively, quarks in a color singlet should have an energy 


( 1 ) 


^ 4 0^5 

E --^ 

3 R 

at very small distances. The figure shows that this distance scale is evidently smaller than 0.1 fm. 

The actual behavior of the adiabatic potentials at small distance is problematic: lattice gauge 
computations diverge at the origin and must be regulated in some way. Furthermore, level crossing 
is expected to occur; namely, it will become energetically favorable for a configuration to emit a 
scalar glueball and convert to a qq color singlet state at some distance. This surface and the level 
crossings are indicated by the dashed line in the figure. Of course the quantum numbers must agree 
for the transition to occur and the mixing matrix element must be large enough to resolve it at the 
temporal extent used in the lattice computation. 

The large distance behavior of the adiabatic surfaces is of interest because it is widely held 
that this behavior is governed by string dynamics [^[^. The simplest geometric string model is 
described by the Nambu-Goto action 




( 2 ) 


where is a massless vector field with two transverse components. In D dimensions this gives rise 
to a spectrum (computed with fixed ends) of 11 (the string tension is denoted (j)[^ 


D-2 


2ttN\ 


1/2 




( 3 ) 


^Note that D = 26 is required for consistent quantization, but this problem is not relevant for large R. 
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Figure 1: Lattice Adiabatic Hybrid Potentials. Figure curtesy of C.J. Morningstar [^. 


The fact that this formula works quite well for the ground state surface for i? > 0.5 fm has been 
traditionally taken as an indication for the robustness of the string description of gluodynamics 


At large distances the Nambu-Goto action predicts splittings between adiabatic surfaces that 
behave as tt/R. Detailed comparison to large distance lattice results by Juge et al reveal that 
the expected universal behavior is achieved, but only for quark separations greater than 2 fm; 
furthermore residual deviations from string behavior indicate the presence of fine structure in the 
interaction, so that an effective string interaction is more appropriate for the description of gluonic 
excitations at large distances [T^ . 


2.2 Gluelumps 


A “gluelump” is a hadronic state comprised of a static source in the octet representation with 
accompanying glue, such that the resulting state is a singlet under gauge transformations. Investi¬ 
gation of gluelumps was motivated by interest in the properties of bound states of massive gluinos 
and gluons 13 . The initial lattice investigation was made by Michael in 1985 in quenched SU(2) 


gauge theory 14 m- 

We note that the gluelump spectrum can be determined up to an infinite self energy due to the 
static adjoint gluon. Thus energy differences can be unambiguously extracted. Absolute energies 
can also be obtained under specified renormalization schemes. Here we present the results of Bali 
and Pineda [^, which were obtained in quenched SU(3) gauge theory. Absolute gluelump energies 
were obtained by employing a renormalon subtraction scheme [Tt] with a matching scale set to 
2.5/ro ~ 1 GeV. Results for the lowest lying eight levels are shown in Table The authors of 
Ref. have examined a large set of quantum numbers and non-octet sources, with results similar 
to those presented here. 
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mass (GeV) 

1+- 

0.87(15) 

1“ 

1.25(16) 

2“ 

1.45(17) 

2+- 

1.86(19) 

3+- 

1.86(18) 

0++ 

1.98(18) 

4“ 

2.13(18) 

1-+ 

2.15(20) 


Table 1: Low Lying Gluelump Spectrum [Tg] . 


HERE 

Because the quark and antiquark can be regarded as merging into a color octet as the separation 
between them lessens, the spectrum of adiabatic hybrid surfaces is related to the gluelump spectrum 
as R 


0 19 


In this limit greater symmetry is obtained because the gluonic spin, Jg^ becomes a 


good quantum number |15[|20| and the gluelump spectrum can be mapped to the adiabatic surfaces 
as shown in Table The match between the adiabatic and gluelump spectra is displayed in Fig. 

where reasonably good agreement is seen. The exception is the which appears headed 

towards the 0^^ gluelump energy, but then deviates downwards. This is perhaps an example of 
the adiabatic surface crossing mentioned above. In fact, the final point on the surface appears 
to have landed on the surface given by R) + mQ++, shown in the figure as a cyan line. 


gluelump 

adiabatic surface quantum numbers 

i+- 


1— 

np, sy 

2— 


2+- 


3+- 

n", A(,i, 

0++ 

Ey' 

4— 

E-', n", A^i, $ 3 , Tg 

1-+ 

V + ' TT'" 


Table 2: Adiabatic Gluon Surface Degeneracies at Small Distance. The 2+ and 3+ identifications 
could be reversed. 

2.3 Heavy Hybrids 

The simplest method to access hybrids with heavy quarks is to leverage the adiabatic gluonic 
potentials of section |2.1| by solving the non-relativistic Schrodinger equation for quark motion on 
the relevant surface. This approach has been adopted by Juge, Kuti, and Morningstar [^, who 
address a non relativistic Schrodinger equation with a centrifugal factor given by 


(LJq) = L(L + 1)-2A2 + (j2), 


( 4 ) 















Figure 2: Adiabatic Hybrid Potentials with Gluelump Spectrum. (Figure reproduced with permis¬ 
sion from reference [Tg]). 


where L is the total angular momentum of the system. The result is a spin-averaged spectrum that 
depends on the adiabatic surface employed as a potential. For the 11+ surface one has (J^) = 0 and, 
since the D+ surface mimics the phenomenological Cornell potential, the usual results for heavy 
quarkonia are recovered. For excited surfaces, H,^ and one can approximate (J^) = 2 since this 
corresponds to the minimum gluonic field angular momentum permitted in these configurations. 

Meson quantum numbers are constructed in terms of the total quark spin the total angular 
momentum L, and the mesonic spin J — L + S. In the leading Born-Oppenheimer approximation 
the eigenvalues L{L -\- 1) and S{S -\- 1) of and S‘^ are good quantum numbers and the parity P 
and charge conjugation C of a meson is given in terms of L, 5, and A by 


P = e(-)^+'^+^ 

(5) 

C = 

(6) 


Here L > A e = ± for and both signs apply for A > 0. Note that r] is the product of gluonic 
charge and parity quantum numbers defined above. Low lying quantum numbers in the ground 
state 11+ state are thus as given in Table 

This is, of course, the usual pattern of the non-relativistic constituent quark model. 

For the phenomenologically relevant H-^^ surface the low lying quantum numbers are (we take 
A = 1, 77 = —1, e = ±) given in Table 

Because the spectrum only depends on the radial quantum number and L at this order, the 
leading Born-Oppenheimer multiplets are (0 ^,1 ), (l+“, (0,1, 2)++), etc, for the 11+ surface 

and 

1--, (0,1,2)-+; 1++, (0,1,2)+- (7) 
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L 


0 

0 

0 -+ 

1 

0 

1 — 

0 

1 

1 +- 

1 

1 

( 0 , 1 , 2 )++ 


Table 3: Meson Quantum Numbers. 


5 

L 


0 

1 

1 —, 1 ++ 

1 

1 

( 0 , 1 , 2 )-+, ( 0 , 1 , 2 )+- 

0 

2 

2 ++, 2 - 

1 

2 

(1,2,3)+-, (1,2,3)-+ 


Table 4: Meson Quantum Numbers. 


for the lowest states on the surface. We remark that these multiplets do not agree with 
those observed in explicit lattice hybrid computations. Rather, the low lying multiplet contains 
— 1 , ( 0 , 1 , 2 ) which corresponds to e = + 1 . 

Fitting the S-wave D+ energy to the spin-averaged 775 and T mass gave a bottom quark mass 
of 4.58 GeV. This was then used to obtain the spectrum shown in Fig. Experimental results 
are given as solid lines in the figure. Notice that the agreement of the conventional spectrum 
with experiment appears to worsen dramatically once past BB (BB*^ S*S*) threshold. Although 
similar large shifts due to light quarks may be present in the hybrid spectrum, there is reason to 
believe that these will be small for the lowest lying states. This is discussed further in section]^ 
Lattice studies with heavy quarks are difficult because heavy quarks tend to lie above the 
ultraviolet cutoff on typical lattices and hence are removed from the dynamics of the theory. One 
way to avoid this problem is to work directly with non-relativistic effective field theory, which 
removes the heavy quarks in an ordered expansion. An early study of this type by the CP-PACS 
collaboration neglected all spin-dependent operators in the effective Lagrangian [^. The authors 
employed a “magnetic” hybrid interpolating fields for a spin singlet BiX and spin triplet 

= ^p^ajBiX states ( 7 /^^ and y are the heavy quark and antiquark creation operators respectively). 
Since spin-dependent operators were not included in the Lagrangian, the result was a degenerate 
multiplet of hybrid states with quantum numbers 1 , (0,1,2) 

The computations yielded a charmonium hybrid multiplet 1.323(13) GeV above the spin av¬ 
eraged charmonium ground state (i.e., 3.069 GeV) and a bottomonium hybrid multiplet 1.542(8) 
GeV above the ground state (at 9.445 GeV). 

A calculation of hybrid masses with large lattices in quenched QCD has been performed by 
Bernard et al 20 . They obtained charmonium hybrid masses of 


1“+ 4.39(8) GeV 0+“ 4.61(11) GeV. ( 8 ) 

The error in these results is statistical only - additional truncation and quenching errors also exist. 

Advances in computer technology, especially leveraging the power of graphical processing units, 
has permitted the recent computation of charmonium states with the full QCD Lagrangian. In 
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Figure 3: Leading Born-Oppenheimer Spectrum for bottomonium (S, P) and Hybrids (H-^^ and ). 
(Figure reproduced with permission from reference [^). 


particular, the Hadron Spectrum Collaboration has performed a large scale unquenched calculation 
that employs a large variational basis, a fine temporal lattice spacing, two light dynamical quarks, 
a strange dynamical quark, and improved lattice actions to obtain a comprehensive charmonium 
spectrum [^. Despite these technical advances, the dynamical quarks are still heavy, yielding 
a pion mass of 396 MeV, and a — rjc splitting of 80(2) MeV - too small compared to the 
experimental value of 117 MeV. 


Despite these shortcomings, the spectrum, shown in Fig. is of relatively high quality. The 
figure shows charmonium masses as boxes (with vertical extent indicating errors) and experimental 
states as lines. One sees good agreement for low lying states, with diminishing accuracy higher in 
the vector channel. Numerical values for the computed masses are presented in Table 


The authors of Ref. 23 


also measured state overlaps with various operators, (M|O|0), as a 
probe of the internal structure of the state \M). Thus, for example, some vectors have large 
overlaps with a quark-antiquark pair in a ^Si state, while others have larger overlap with ^Di 
operators. These overlaps only provide qualitative indications of state configurations because they 
are scale-dependent and comparison to continuum matrix elements can be confounded by operator 
mixing. Nevertheless, comparison with conventional quarkonia structure as determined by the 
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Figure 4: A lattice QCD calculation of charmonium States. (Figure reproduced with permission 
from reference [^). 


non-relativistic quark model [24[j^ reveals that the method is reliable. 

This method can be used to determine states with large overlaps with operators with gluonic 
content. The resulting states are indicated with red and blue boxes in Fig. As can be seen, an 
approximate multiplet forms with quantum numbers 1 and (0,1, 2) This structure can easily 
be obtained if the effective gluonic degrees of freedom have = 1^ . Combining this with a 

qq pair in a state yields the 1 component of the hybrid multiplet, while combining it with 
qq in ^Si yields the remaining members. Higher quantum numbers obtained in this manner are 
listed in Table. Note that the higher hybrid multiplet in Fig. [^contains — (0,1, 2,)++ and 
(0,1, 2, 3)^ states, which maps very well to the expected P-wave multiplet. 

2.4 Light Hybrids 

Computations with light quarks are more difficult than those with moderate mass quarks because 
the Dirac matrix that must be inverted becomes rapidly ill-conditioned and larger (a lighter pion 
requires a larger lattice to maintain constant physics). Furthermore, light quarks mean that decays 
to multiple pion states are possible, greatly complicating the extraction of observables. 

The earliest computations of light hybrid masses were made in the quenched approximation 
which effectively ignores internal qq loops. These calculations all predicted that the 1 ^ nonet of 
hybrid mesons was the lightest, with a masses in the 1.8 to 2.1 GeV mass region [27[j^ . 

A summary of 1 ^ hybrid mass computations is shown in Fig. Open symbols represent 
quenched computations, while filled symbols are unquenched. General agreement in the data is 
evident; and a naive extrapolation to the physical pion mass gives an unquenched 1 ^ mass of 
approximately 1.6 GeV. 

Finally, we present the results of a recent exhaustive computation of the light meson spectrum 
by the Hadron Spectrum Collaboration [^. This computation was made with a large operator 
basis, on lattice of size 16^ x 128 up to 24^ x 128 lattice, a temporal lattice spacing around 0.034 fm 
and a spatial lattice spacing of approximately 0.12 fm, and pion masses of 702, 524, and 391 MeV. 
Hairpin (disconnected) diagrams were included with the aid of the “distillation” method (^. The 
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(M - / MeV 

0-+ 

0 663(3) 1143(13) 1211(13) 

1“ 

80.2(1) 698(6) 840(3) 1154(28) 1301(14) 1339(38) 

2-+ 

860(3) 1334(17) 1350(17) 

2“ 

859(5) 1333(18) 

3“ 

867(3) 1269(26) 1392(12) 

4“+ 

1444(10) 

4“ 

1427(9) 

0++ 

461.6(7) 972(9) 1361(46) 1488(30) 

1+- 

534(1) 1006(9) 1360(38) 1462(51) 1493(19) 1513(39) 

1++ 

521.6(9) 1002(10) 1415(14) 1484(48) 

2++ 

554(1) 1041(12) 1112(8) 1508(21) 

3+- 

1142(6) 1564(22) 

3++ 

1130(9) 

4++ 

1129(9) 

1-+ 

1233(16) 

0+- 

1402(9) 

2+- 

1411(40) 1525(18) 


Table 5: Charmonium Spectrum 
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1+- 

0 

0 

1— 

1+- 

0 

1 

(0,1,2)-+ 

1+- 

1 

0 

(0,1,2)++ 

1+- 

1 

1 

(0,1,1,1,2,2,3)+- 


Table 6: Hybrid Multiplets 


computation did not include glueball or multihadron operators, and thus the extracted “masses” 
are only approximations to the resonance mass parameters. 

Results for the isovector and isoscalar spectra are shown in Fig. for pion masses of 392 
MeV. Notice that mixing between light and strange quarks is represented (in green and black) in 
the figure. States outlined in orange have large overlap with gluonic operators. Notice that the 
quantum number-exotic states (to the right) are all predicted to be approximately ideally mixed. 

The quark mass-dependence of the lightest hybrid multiplet is of phenomenological interest. 
This has been evaluated by the Hadron Spectrum Collaboration and is discussed by Dudek in 
Ref. (^. Fig. 0 shows the multiplet for four different pion masses. It is apparent that the P-wave 
0^ and 2^ are approximately independent of quark mass, implying that a short range poten¬ 
tial dominates the effective hybrid spin-dependent potential. Furthermore, the 1 state is also 
largely spin-independent, which implies that quark spin-triplets are required in the spin-dependent 
potential. Finally, it appears that the spin-triplet J ^ multiplet slowly splits as the quark mass 
is reduced, with the J = 0 component decreasing slowly, the J — 1 decreasing more rapidly, and 
the J — 2 increasing slowly. Developing a phenomenological model for these observations is an 
interesting task. 
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Figure 5: The mass of the — 1 ^ exotic hybrid as a function of the pion mass from lattice 
calculations. The open (cyan) symbols correspond to quenched calculations, while the solid (red 
and blue) symbols are dynamic (unquenched) calculations: open (cyan) star open (cyan) 
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2.5 Hybrid Baryons 

All quantum numbers are available to baryons, hence ‘exotic’ quantum number baryons do not 
exist (perhaps explaining the relative lack of interest in these states). Furthermore, conventional 
and hybrid baryons will mix to form the physical spectrum, which can seriously affect the ability 
to compute the properties of these states and to discover them experimentally. 

We are aware of only one lattice computation of the hybrid baryon spectrum, which was carried 
out by the JLab group [^. The authors considered the spectrum of nucleons and deltas at several 
quark masses and found a set of positive parity hybrid baryons with quantum numbers 2 [A^x/ 2+]5 
2 [A/' 3 / 2 +], A^ 5 / 2+5 ^ 1 / 2+5 2 ind A 3 / 2 + above the first band of conventional excited positive parity 
baryons. Results are shown in Fig. These have been obtained on a 16^ x 128 lattice with 2+1 
dynamical quarks and a pion mass of 396 MeV. 

The low lying hybrid baryons appear as a band with positive parity that lie above the first- 
excited positive parity conventional states. This pattern of states is compatible with a color octet 
gluonic excitation with = 1^, in keeping with hybrid mesons. Furthermore, the hybrid excitation 
scale is approximately 1.3 GeV for both mesons and baryons, which indicates a common mechanism 
in both of these systems. It is notable that hybrid states appear to mix weakly with conventional 
states, which encourages the idea that they might be experimentally detectable. 


2.6 Transitions 

Relatively little lattice work has been done on hybrid transitions. Indeed, the difficult task of 
evaluating two-point correlation functions such as ( 02 (^ 2)01 (+)) and extracting physical states is 
replaced by the daunting task of evaluating three point functions { 02 {t 2 )ji^{ti) 0 i{ti))^ extracting 
physical states and normalizing correctly. 
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Figure 6: Isoscalar and Isovector Hybrid Spectrum. States outline in orange have large gluonic 
content, (Figure reproduced with permission from reference [38]). 
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Figure 7: The Lightest Hybrid Multiplet as a function of Quark Mass. (Figure reproduced with 
permission from reference [^). 


2.6.1 Vector Hybrid Mixing 

Computing an overlap matrix element is a somewhat simpler task than measuring a transition 
matrix element. For example, the MILC collaboration has examined the effects of moving beyond 
the static limit in the heavy quark expansion of quenched NRQCD by evaluating the matrix element 
{QQ\0\QQg) [^. The operator is the leading order correction in the NRQCD Lagrangian: 


0 = -cg 


a-B 
2M ' 


(9) 


where M is the heavy quark mass, B is the chromomagnetic field, and c is a coupling that is unity 
at tree level. 

The resulting matrix elements can be written in terms of admixture fractions as 


|(T|QQ5(1--))|2«0.4%, 

(10) 

I(%IQ05(0-+))|2«i%, 

(11) 
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Figure 8: Nucleon and Delta States. Grey boxes are conventional baryons; blue boxes are hybrid 
baryons. (Figure reproduced with permission from reference [^). 


\{J/^P\QQg{l—))f^2.3%, 

(12) 

\{vc\QQg{o~^))\‘^ ^ 6%. 

(13) 


The authors caution that the results for charmonium are less reliable than those for bottomonium 
since higher order operators can contribute substantially. Furthermore, unquenching effects can be 
significant. Neglecting these issues permits a simple estimate for the rate of charmonium vector 
hybrid production in electron-positron collisions: 

V{ccg e+e“) 0.023(2) • V{J/^ e+e“) ^ 0.12(1) keV. (14) 

2.6.2 Heavy Hybrid Hadronic Transitions 

The first lattice computation of a hadronic transition was made by the UKQCD collaboration for 
the case of heavy hybrids [^. The static quark limit imposes important constraints on the decay 
process since the quark-antiquark configuration must remain unchanged. The authors focus on the 
decay of the exotic 1 ^ state and determine that decay into S-wave mesons is forbidden (since 
production of the light quark pair in a spin triplet is forbidden by r] while a spin singlet is forbidden 
by A). Furthermore, decay to an S-wave {Qq) + P-wave (qQ) configuration is forbidden because 
the P-wave excitation energy is typically greater than the hybrid excitation energy. Thus the only 
allowed transition in the heavy quark limit is a hidden flavor process, QQg QQ + S', where S 
is a light flavor singlet meson. This amounts to a de-excitation of the excited gluonic string by 
emission of a light quark-antiquark pair. 

The authors computed two such transitions using unquenched QCD with light quark masses 
near the strange quark mass. Heavy valence quarks were fixed in place and integrated out, so that 
they did not appear explicitly in the computation. Results (applied to bottomonium) were 

bbg{l ^) ^ 77 ^ 77 ( 55 ) ^ 1 MeV (15) 
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and 


bbg{l ^ Xb^{ss) ^ 60 MeV. 


(16) 


2.6.3 Light Hybrid Hadronic Transitions 

The decays of light hybrids are much less constrained than those with heavy quarks, thus lat¬ 
tice computations become even more important to guiding experiment and phenomenology. The 
UKQCD collaboration has examined the decay of a light 1 ^ exotic (termed the tti) with two 
dynamical quarks [^. They determined the tti mass to be 2.2(2) GeV and obtained effective 
couplings for two decay modes as follows 

r(7ri ^ biTv)/k = 0.66 ± 0.20 (17) 

r(7ri ^/i7r)/fc = 0.15 ±0.10, (18) 

where k is the relative momentum in the final state. Assuming that these results determine a 
constant effective coupling then permitted the authors to obtain partial widths of 


r(7ri ^ biTv) = 
r(7ri ^ /i7r) = 


400 ± 120 MeV 
90 ± 60 MeV. 


(19) 

( 20 ) 


As a check of their procedure, they also carried out a similar calculation for bi cjtt, obtaining 
T/k ^ 0.8, which leads to T{bi uott) ^ 220 MeV. This is about a factor of 1.6 larger than the 
experimental width. Finally, the authors noted that the large hybrid widths arose chieffy due to the 
large available phase space in the decay. We remark that in general this expectation is false because 
hadronic form factors cause effective couplings to diminish rapidly with increasing momentum. 

The latter point was subsequently pursued by Burns and Close, who compared the lattice results 
for the transition element near threshold to those of the ffux tube model (see section 4.1) 45 . The 


two approaches were found to be in rough agreement near threshold, see Table Thus, if the ffux 
tube model can be trusted to extrapolate to the physical recoil momentum, one obtains substantially 
reduced partial widths of T^tti bin) ^ 80 MeV and T{7Ti fin) ^ 25 MeV. The lattice results 
also suggest that the light quark creation vertex has spin triplet quantum numbers. 



bi ujn (GeV ) 

TTi ^ bin (GeV ^/^) 

TTi ^ fin (GeV ^/^) 

Lattice (UKQCD) 

2.3(1) 

2.9(4) 

1.5(4) 

Lattice (CP-PACS) 

3.4(2) 

2.9(3) 

1.1(4) 

ffux tube (a) 

2.7 

2.9 

1.4 

ffux tube (b) 

3.3 

3.9 

1.9 


Table 7: Comparison of lattice and ffux tube model transition amplitudes. UKQCD and CP-PACS 
refer to different gauge configuration ensembles. 


2.6.4 Hybrid Radiative Transitions 

More recently, the Hadron Spectrum Collaboration has computed charmonium hybrid radiative 
transitions [^. The calculation was made with a large operator basis in the quenched approxi¬ 
mation. The renormalization constant required to compare the lattice matrix elements to physical 
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ones was determined nonperturbatively by conserving charge at zero recoil. The resulting widths 
are presented in Table where one sees quite good agreement with experiment, where available. 
Notice that the process ccg{l ^ is a magnetic dipole transition. With conventional 

charmonia, these require a spin flip and therefore are suppressed for heavy quarks. In this case, 
however, the extra gluonic degrees of freedom can permit the transition, and hence it can be large. 


transition 

rlattice (keV) 

Texpt (keV) 

XcO ^ J/'fpl 

199(6) 

131(14) 

'>P' Xc07 

26(11) 

30(2) 

V’" ^ Xco7 

265(66) 

199(26) 

ccg{l ) Xc07 

< 20 


J/i> Vcl 

2.51(8) 

1.85(29) 

V’' rjc'y 

0.4(8) 

0.95 - 1.37 

V Vcl 

10(11) 


ccg{i ) ^ Vci 

42(18) 


ccg{l ^ 

115(16) 



Table 8: Quenched Lattice Charmonia Radiative Decays [4^ . 


3 Hybrid Models 

The construction of a reliable model of hybrid meson structure and dynamics is important for 
the interpretation of experimental results. At the simplest level, this is because it is expensive to 
compute large numbers of experimentally relevant quantities on the lattice. It is also likely that 
the computation of complicated amplitudes involving hybrids will remain out of reach of lattice 
methods for a long time. 

In this section we discuss the salient features of bag models, string models, constituent glue 
models, and attempted computations with QCD sum rules and the Schwinger-Dyson formalism. 
The prime difference between models is the assumed form that the gluonic degrees of freedom take 
on: broadly, these are quasiparticle or collective in nature. 


3.1 Bag Models 


A detailed phenomenology of hybrids was first developed with bag models, both in the MIT bag 
model (47-52 , and in the “Budapest variant” 53 . The idea was to place a gluonic field in a 


vacuum Tubble’ with appropriate boundary conditions. These imposed either transverse magnetic 
(1 , TM) or transverse electric (1^ , TE) solutions, with the TE modes being the lightest. The 

predicted result was four low lying nonets with, in order, quantum numbers 0 1 1 , and 

2 ^ in the mass range 1.2 - 2.5 GeV. 

The Budapest variant was designed to deal with heavy quarks for which a spherical bag was 
unrealistic (as it would not naturally yield linear confinement). Thus the bag was allowed to deform 
due to the presence of fixed heavy quark and antiquark sources. The resulting adiabatic energy 
surface was used in a two-body Schrodinger equation to give mass estimates for hybrids. Masses 
found for the lightest hybrids were 3.9 GeV for cc and 10.45 GeV for bb. Readers interested 
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in examining this model closely should be aware that the formulas for quantum numbers reported 
therein are incorrect IM. 


The bag model has also been applied to hybrid baryons . In this case a TE gluonic field was 
combined with three quarks in an overall color octet to produce a set of hybrid states: ‘ 2 [Ni/ 2 ^]^ 
2[A^3/2+], A^ 3 / 2+5 A^ 5 / 2+5 ^ 1 / 2+5 ^ind A 3 / 2 +- Notice that these quantum numbers agree with those 
reported for the lattice in section [23 


The overall mass scale of these states is not firmly established; for example in Ref. the scale 
was set assuming that the 7r(1300) was a hybrid meson, yielding hybrid baryons near 1.6 GeV. 
These days the 7r(1800) is a more plausible hybrid candidate, which places the low lying bag model 
hybrid baryons above 2 GeV in mass. 


We remark that although bag models provide a compelling qualitative picture of confined quarks 
and gluons, they suffer from a number of ambiguities which makes their application less than ideal. 
Among these are difficulties in determining the response of the bag boundary when quarks and 
gluons are present, ambiguities in gluon self-energies, and the existence of spurious degrees of 
freedom associated with the center of mass. 


A more subtle problem concerns the presence of both a bag pressure and valence gluons. This 
issue is illustrated with a statement by Jaffe and Johnson, who discuss exotic quantum numbers in 
the MIT model and note that they can occur due to “type II exotic quark states”, which are qq or 
qqq states with quantum numbers that cannot be obtained in the “conventional” quark model 56 . 
They contrast these with “mixed states of quarks and glue” and go on to explain that it is the bag 
surface that permits the extra quantum numbers. 


The problem is that the degrees of freedom in QCD are quarks and gluons, and it is natural 
to associate the bag pressure with gluonic properties of the medium. Thus the extra quantum 
numbers are due to glue, and “type 11” and states are really the same. In short, an intrinsic 

ambiguity exists between bag and gluonic degrees of freedom. 


More modern applications of bag models have focussed on reproducing lattice results. For 
example, Juge, Kuti, and Morningstar have resurrected the Budapest variant and compared its 
adiabatic surfaces to those of the lattice [^. The results, some of which are given in Fig. show 
quite good agreement between the approaches. 


Finally, an attempt at reproducing the gluelump spectrum in the bag model has been made 
by Karl and Baton [^. The calculation draws on the mapping between adiabatic surfaces at zero 
quark separation and the gluelump spectrum. Computed single gluon modes in a spherical bag do 
not agree with the lattice gluelump ordering; including corrections due to the Coulomb interaction 
and adjusting the strong coupling gave levels as reported in Table For comparison we display 
the lattice results of table [l] with energies shifted to match 1.43 GeV for the 1^ state. As can be 
seen, there is broad agreement, while detailed agreement cannot be claimed. 


J^: 

1+ 

2- 

1“ 

3+ 

2+ 

mass (GeV): 

1.43 

1.97 

1.98 

2.44 

2.64 

lattice (GeV): 

- 

2.01 

1.81 

2.42 

2.42 


Table 9: Gluelump Bag Model Predictions and Shifted Lattice Spectrum. 
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Figure 9: Lattice and Bag Model Adiabatic Potentials. (Figure reproduced with permission from 
reference [^). 


3.2 Relativistic String Models 

Models of soft glue as collective, ‘stringy’, degrees of freedom have a long history, dating from early 
ideas contained in dual resonance models. The formal development of string theory veered off into 
mathematical details once it was realized that substantial formalism was required for consistent 
quantization. However, variant string models were presented as explicit models of mesons. Thus, 
for example, Andreo and Rohrlich noted that dual strings and gauge theories are related, and thus, 
“one is ... led to believe that in a certain approximation a meson can be regarded as two point-like 
quarks confined to each other by a string” [^. Once this idea is in place it is natural to examine the 
adiabatic potentials of the model and note that “daughter trajectories” correspond to additional 
mesonic states. 

To our knowledge, the first to draw this connection were Giles and Tye in 1976 [^. These 
authors coupled quarks to a relativistic two-dimensional sheet in what they called the “Quark 
Confining String Model” and noted that “The presence of vibrational levels gives ... extra states 
in quantum mechanics ... that are absent in the charmonium model”. 

The quark-string system was assumed to obey an action 

5 = J (-det(5))i/2 + ... (21) 

where r is a collection of tangent vectors to the string surface. Vibrational energies were crudely 
estimated and a Born-Oppenheimer Schrodinger equation was solved for the hadronic states: 
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with 


H = - 


1 

m dB? 


+ 2m + Vn{R) + 


mB? 


( 22 ) 


K = fei? 1 - 


2n7r 


2n7r + A:((i? - 2d)2 + 4^2) 


- 1/2 


(23) 


Here is a correction due to the finite quark mass. Vibrational levels (hybrids) in the charmonium 
sector were predicted to lie at 3.96, 4.21, 4.41, 4.45, and 4.46 GeV; although no quantum numbers 
(other than the angular momentum) were given. 

Similar models have been developed by a variety of groups: Ref. extended the work of Giles 
and Tye and predicted a vector bottomonium hybrid with mass mT + 990 MeV and a photocoupling 
of V(hhg{l ) ^ e+e“) = 0.20(15) keV. An early application of string models to baryons is 
contained in Ref. where adiabatic potentials relevant to the three quark system are developed. 
Another application has been to the adiabatic surfaces of section 2.1 The authors of Ref. used 


a simple quantized string model with fixed ends to obtain an expression for the adiabatic surfaces 


Vn{R) = y(ai?)2 + 2na{N + c) + C (24) 

where N labels the string excitation and is given in terms of phonon occupation numbers as 


V = 


oo 

E 

m=l 


m{n, 


m “1“ ^m) 


(25) 


The authors added —4/3 • ag/B to this expression, set c = 0, and fit to the adiabatic 11+ surface to 
obtain C and a. The resulting excited surfaces agreed reasonably well with the lattice computations. 

The extension of effective string models to 1/M^ corrections to the ground state adiabatic 
potential has been explored by Brambilla et al [^. Explicit expressions for a variety of corrections 
were obtained, but a detailed comparison to quarkonium phenomenology was not made. 

Recently, the AdS/QCD hypothesis of Maldacena has been used to extract the ground state 
and the 


66 


E+ surface 
with lattice results is quite good. 


hybrid surface 67 . Although this effort is in its infancy, the agreement 


3.3 Flux Tube Model 
3.3.1 Model Development 

The flux tube model was developed by Isgur and Patou as a response to bag models. In particular, 
Isgur and Patou objected to the use of perturbative gluons in constructing a hadron, and argued 
that nonperturbative degrees of freedom, such as in a string model, should be used 68 . Noting 


that the QCD lattice Hamiltonian also is not written in terms of perturbative gluons, Isgur and 


Patou built a simple model by truncating the Hamiltonian with a series of approximations [69,70 


In particular, the lattice degrees of freedom are quarks fields on lattice sites and gluonic dink 
variables’ Ui = exp(—ia^A^(x)) where £ represents the link (x,/i). In the strong coupling limit the 
Hamiltonian is given by 
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( 26 ) 


HscQCD — 2^ ^ ^ ^ai ^ ^ ^ 
i n 

where g is the strong coupling, a is the lattice spacing, and n is a lattice site. The velocity 
variables Ui have been replaced by electric field operators E^. Gauge invariant pure glue states are 
formed by closed (possibly multiply connected) loops of link operators. The commutation relation 
U] = T^U then implies that the energy of these states is simply the sum of the quadratic color 
charges of each link: 


E\oop — 2^ Ci (27) 

^Gloop 

where = 4/3 for a field in the 3 or 3 representations, 10/3 for 6 or 6, etc. The presence of quarks 
permits gauge invariant states with open flux strings which terminate on quark color sources or 
sinks. Perturbations to these states are provided by subleading quark hopping and magnetic terms. 
The former allow flux tube breaking via quark pair production or quark motion. The latter can 
change link color representations, cause link hopping, or change loop topology. 

Isgur and Patou simplified the dynamics by (i) assuming an adiabatic separation of quark and 
gluon degrees of freedom (ii) neglecting ‘topological mixing’ such as loop breaking or loop Euler 
number changing transitions (iii) working in the non-relativistic limit. The model is meant to be 
applied to the ‘intermediate regime’ where 1/a ^ Vb. They then modelled link variable dynamics 
in terms of spinless colorless particles (‘beads’) of mass ba where b is the string tension in the static 
quark potential. Finally these particles are assumed to interact via a linear potential and perform 
small oscillations about their resting positions. The result is a simple discrete string model for glue 
described by the Hamiltonian: 


H = bQR + J2 

n 


pI .bo 2 

^ + ^(yn-yn+i) 


(28) 


where yn is the transverse displacement of the nth of N string masses, pn is its momentum, bo 
is a bare string tension, and R = {N + l)a is the separation between the static quarks. This 
Hamiltonian may be diagonalized in the usual way yielding 


HfTM - boR + (:^-R ~ ~~ + • • •) + (29) 

where creates a phonon in the nth mode with polarization A. Notice that the string tension has 
been renormalized by the first term in the brackets. The last term in the brackets is the Liischer 
term of string phenomenology [^. The mode energies are given by uon — 2/asin[7rn/2(A^ + 1)]. 
Thus uji ^ 7r/i? as ^ (X) is the splitting between the ground state Coulomb+linear potential 
and the first gluonic excitation surface at long range. The energy of a given phonon state is 
approximately 


E — Eo + ^ ^ 

with N given by Eq. 


(30) 
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Hybrid mesons are constructed by specifying the gluonic states via phonon operators and com¬ 
bining these with quark operators with a Wigner rotation matrix: 


|LML;ss;A,{n„^+,n^_}) (X f (31) 

^ ’ ’ m 

The projection of the total angular momentum on the qq axis is denoted by A = ~ 

The parity and charge parity of these states are given by 

P\LML]SMs]K{nm+,n^-]) = -A, {^m-, n^+}), (32) 

C\LML\SMs\K{nm+,n^-]) = (-)^+'^+^+^|LMi; 5Ms;-A, {n^-, n^+}). (33) 

States of good parity are thus formed as 


\LMl] SMs\ C; a, {nm+,n^-}) = {\LMl] SMs] A, {nm+, n^_})+C|I/ML; SMs', -A, {n^_, n^+})^. 

(34) 

Possible single phonon (m = 1) mesons are listed in table 10, where underlined quantum 
numbers represent quantum number exotic hybrids. 


c 

L 

5 


+ 

1 

0 

1++ 

+ 

1 

1 

(2,1,0)+- 

+ 

2 

0 

2— 

+ 

2 

1 

(3,2,1)-+ 

- 

1 

0 

1— 

- 

1 

1 

(2,1,0)-+ 

- 

2 

0 

2++ 

- 

2 

1 

(3,2,1)+- 


Table 10: Flux Tube Model Single Phonon Mesons 


Isgur and Paton obtained hybrid meson masses by solving a model Hamiltonian of quark motion 
on the single-phonon excited surface: 


HiP = 


y ^ A(L + l)-A^ _ _ e-f^bR 

2ndB? 2nB? 3R 


(35) 


The interaction term incorporates several important additional assumptions. Namely the tt/R 
phonon splitting is softened at short range. The parameter / which appears in the softening 
function was estimated to be roughly unity. Furthermore, it was assumed that the attractive 
Coulomb (1/i?) potential remains valid for hybrid mesons. 

Finally the quark angular momentum operator is now complicated by the presence of glu- 
onic/string degrees of freedom. One may write 


Lqq = L - - Ls^ 


(36) 
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where L (Ls) is the total (string) angular momentum. Note that Ls^ 

Using = AR and neglecting surface mixing yields the centrifugal term of Eq 
this is not the same as used by Juge et al). 


mixes adiabatic surfaces. 

(notice that 


35 


The hybrid masses obtained by solving Eq. 35 are labelled Ejp in table 11 Isgur and Eaton 


also estimated the effects of adiabatic surface mixing and used these as their final mass estimates 
(labelled Ejp). The column labelled KW lists hybrid masses obtained when the adiabatic surfaces 
of section |2.1| are used along with the centrifugal term of Eq. [3 . 


ffavor 

Eip (GeV) 

E'lp (GeV) 

Ekw (GeV) 

1=1 

1.67 

1.9 

1.85 

1=0 

1.67 

1.9 

1.85 

ss 

1.91 

2.1 

2.07 

cc 

4.19 

4.3 

4.34 

bb 

10.79 

10.8 

10.85 


Table 11: Hybrid Mass Predictions 


3.3.2 Further Developments 

In developing the ffux tube model, Isgur and Paton assumed adiabatic separation of bead and quark 
motion and that the string executes small oscillations. Both of these assumptions were tested 
by numerically solving a bead-quark-antiquark system with a quantum Monte Carlo algorithm 
[ 73 ] . The results indicate that the small oscillation approximation is accurate for long strings 
but overestimates gluonic energies by an increasing amount as the inter quark distance shrinks. 
Typical energy differences are order 100 MeV at 1 fm. Similarly, the adiabatic approximation 
underestimates true energies by roughly 100 MeV, with slow improvement as the quarks get very 
massive. It thus appears that these approximation errors tend to cancel each other, leaving the 
mass estimates of Isgur and Paton largely intact. 

The effects of adiabatic surface mixing were examined by Merlin and Paton by considering the 
full quark-bead system [^. Although the effects can be quite complicated, with mixing between 
all surfaces possible, they found that the majority of the effects can be absorbed in a redefinition of 
the hybrid potential by including the rigid body moment of inertia of the string in the centrifugal 
term: 


1 1 
2/ii?2 ^ 2iiR^ + ’ 


(37) 


and with a more important effect that increases the strength of the n/r splitting as r becomes larger 
than niqlb. An explicit computation revealed found mass shifts of order -100 MeV for conventional 
S-wave light quark mesons and -h200 MeV for light quark hybrids. 

Merlin and Paton also examined spin orbit forces in the context of the ffux tube mode [7^ . The 
idea was to map the operators of the leading spin orbit term in the heavy quark expansion of QCD, 
namely Vso = ^/2m • cr • S, onto phonon degrees of freedom. Explicit computations revealed that 
spin orbit splittings due to Vso are small and that the majority of the splittings arise from Thomas 
precession, Vth — l/4(^g x '^q) • cr* This was modelled by including the effects of phonons on the 
quark coordinate. The resulting mass splittings for light hybrids are listed in table 12, One sees 


24 













that the lowest member of the octet of light hybrids is predicted to be the 2+ while the heaviest 
is the 0^ . These results do not agree with those of lattice gauge theory. 



2+- 2“+ 1“+ 0“+ 

1+- 

0+- 

1++ 1— 

5M (MeV) 

-140 -20 20 40 

140 

280 

0 0 


Table 12: Spin Orbit Hybrid Mass Splittings. 


Isgur pointed out that the energy carried by the flux tube will change several features of the 
naive quark model (see also [^). For example, zero point oscillation of the flux tube about 
the interquark axis will induce transverse fluctuations in the quark positions, something which is 
not present when the flux tube is treated as a potential. The additional fluctuations have the effect 
of increasing the charge radius of a heavy-light meson (qQ): 



rrin 


rriq + rriQ 


2 

+ 


2b 

TT^rn^ 


C(3) 


(r^) 


(38) 


where the second term in the bracket is the new contribution. Isgur estimated this to give rise to 
a 50% increase in charge radii of light quark hadrons. 


3.4 Hybrid Baryons in the Flux Tube Model 


Capstick and Page have made a detailed study of baryon flux tube dynamics 77 , 78 . This is 


a technically challenging problem due to the multitude of vibrational and rotational modes that 
are available to a Y-shaped string system. However, they have found that the problem simplifies 
considerably because the string junction decouples to good accuracy from the rest of the bead 
motion. Thus a hybrid baryon can be approximated by three quarks coupled via linear potentials 
to a massive “junction bead”. The dynamics of this system are completely specified by the flux 
tube model and variational calculations indicate that the lowest lying hybrids are and 

states at approximately 1870 MeV. Notice that the lattice confirms these lightest states, but that 
the lattice mass is somewhat higher at 2550 MeV (section [2^ 


We remark that lattice investigations of the static baryon interaction have been carried out 79 
The chief point of interest is whether the expected flux tubes form into a W’ shape or a ‘A’ 
shape. This may be addressed by carefully examining the baryonic energy in a variety of quark 
configurations. Current results are mixed, with some groups claiming support for the two-body 
hypothesis and some for the three-body hypothesis [^. Finally, a strong operator dependence 
in the flux tube profiles has been observed [^, which clearly needs to be settled before definitive 
conclusions can be reached. 


3.5 Constituent Glue Models 


Constituent gluon models for hybrids were introduced by Horn and Mandula 83 and were subse¬ 
quently developed by Iddir et al and Ishida et al [84|-[8^. Since these models assume a diagonal 


gluon angular momentum ig their predictions for quantum numbers differ somewhat from other 
models. 

Horn and Mandula assumed a massless = 1“ gluon interacting via linear potentials with 
quarks. The short range repulsive interquark potential was absorbed into a constant with the 
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argument that repulsive forces do not produce vacuum polarization instabilities. Their model thus 
took the form: 


H = 2mc H- — -\- Pg -\- G{\rq 

rUr. 


rgl + \rq-rg\) + C 


(39) 


where the string tension was taken to be G = 0.30 GeV^. 

For the lightest hybrid states (with ig = 0) Horn and Mandula predicted nonexotic quantum 
numbers equivalent to P-wave qq states, since the gluon has = 1“. Exotic quantum numbers 
including 1 ^ are predicted in the higher-lying {iqq^ig) = (1,0) and (0,1) multiplets. Detailed 
spectroscopic predictions for hybrids have not been published using these constituent gluon models, 
and the estimated masses are assigned large uncertainties. A typical result, due to Ishida et a/., is 
1.3-1.8 GeV for light nonexotic hybrids and 1.8-2.2 GeV for light exotics. 

More recently, the idea of constituent glue models has been revived in the context of QCD in 
Coulomb gauge [^. The model takes the Hamiltonian of QCD in Coulomb gauge as its starting 
point. A nontrivial mean field vacuum Ansatz is used to construct self-consistent constituent 
quarks and gluons, and these are used to construct hadrons in the Tamm-Dancoff approximation. 
In this picture the gluon remains transverse, but has a dynamically generated mass. Interactions 
were truncated at the two-body level. An eigenvalue equation was derived for the case of static 
qua rks and the adiabatic energy surfaces were obtained. Although they mimic those of section 
2.1 reasonably well, the level ordering was incorrect. This problem can be traced to the = 1“ 
assignment for the constituent gluon, and hence applies to all models that make this assumption. 

It is also possible to reproduce the gluelump spectrum with this method by setting the interquark 
distance to zero. In this limit the gluonic angular momentum becomes a good quantum number 
and T-parity is simply given by P(—)-^^. Thus degeneracies appear in the gluelump spectrum, 
precisely as shown in section [2^ (although the level ordering remains incorrect). With this method 
the splitting between the lowest two levels was predicted to be 500 MeV, which should be compared 
to the 360 MeV measured in lattice gauge theory (see table [^. 

Despite these failings, it would be premature to dismiss a constituent gluon picture of hybrids 
because important three body interactions exist in the system and these cannot be neglected in a 
nonperturbative context. Fig. [T0| shows some of these interactions; the first two are self-energies, 
the second pair were considered in older constituent models, and (e) and (f) represent important 
additional interactions. In fact, Szczepaniak and Krusinski have shown that the three-body in¬ 
teractions are sufficient to invert the naive parity ordering in the gluelump and adiabatic surface 
Detailed comparison to lattice results show reasonable agreement, although the model 


88 


spectra 

levels tend to be several hundred MeV high. 

A follow-up paper generalized the Coulomb gauge model by fitting the gluonic quasiparticle 
dispersion relation to the gluelump spectrum (rather than obtaining it by minimizing the vacuum 


energy) 1^. The result and lattice data are shown in table 13 



1+ 

1“ 

2- 

2+ 

3+ 

0+ 

4- 

model (GeV) 
lattice (GeV) 

0.87 

0.87(15) 

1.27 

1.25(16) 

1.47 

1.45(17) 

1.59 

1.86(19) 

1.99 

1.86(18) 

1.98(18) 

2.49 

2.13(18) 


Table 13: Lattice and Constituent Gluon Model Gluelump Spectra. 

The modified Coulomb gauge model was also employed in a computation of the heavy hybrid 
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Figure 10: Interactions in a Constituent Hybrid Model. 


spectrum [^. General agreement with lattice results were found, although the predictions appear 
to be systematically high. The 1 ^ charmonium hybrid was predicted to lie at 4.47 GeV. 


3.6 Sum Rule and Schwinger-Dyson Calculations 

QCD sum rules have been used since their inception for the study of nonperturbative properties 


of QCD 91 . Their application to hybrid mesons dates to the beginnings of the field. A collection 
of results for the light quantum number exotic hybrid are presented in table Evidently mass 


estimates have been climbing over the years, with more recent results tending to agree with lattice 
computations. 


1 + mass (GeV) 

be: 

year 

1.3 

92 

1983 

^ 1 

93 

1986 

1.7(1) 

94 

1987 

1.65(5) 

95 

2000 

1.81(6) 

96 

2009 


Table 14: Sum Rule Hybrid Mass Estimates 


Sum rules have also been applied to hybrid baryons in Ref. 97 . The authors considered the 


possibility of a hybrid baryon as the first excited state above the nucleon with A^i/ 2 + quantum 
numbers and concluded that such a state should have a mass near 1.5 GeV. This is 1 GeV lighter 
than current lattice results. 

Finally, exploratory computations of hybrid masses have been carried out with the Schwinger- 
Dyson and Bethe-Salpeter formalism 98 , [^. Although the authors solve a two-body Bethe- 


Salpeter equation, exotic quantum numbers are accessible because the formalism is relativistic and 
the Bethe-Salpeter kernel can contain gluonic degrees of freedom. Of course the dynamics and 
eigenvalues of hybrids then depends crucially on the assumed form of the scattering kernel. In this 
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case, a separable kernel Ansatz of the form 

Kip,q) = +p-qF{p^)F{q‘^)]r (40) 

was used. The functions F and G are related to the full quark propagator via the Schwinger-Dyson 
equations for the gluon propagator. Results for a selection of states are presented in table along 
with possible resonance assignments suggested by the authors. The authors point out that the 1 ^ 
states carry negative norm, which is normally regarded as a sign of a spurious state although 

they argue that this is not enough to dismiss the states as physical. 



mass (GeV) 

assignment 

0++ 

749 

(j(540), ao(980) 

0+- 

1082 

exotic 

0“ 

1319 

exotic 

1“ 

730 

p(770) 

1+- 

1244 

/ii(1170) 

1++ 

1337 

ai(1260) 

1-+ 

1439, 1487 

7ri(1400), 7ri(1600) 


Table 15: Bethe-Salpeter Hybrid Masses 


3.7 Model Comparison 

It is instructive to confront the models of the past section with the lattice data for the gluelump 
spectrum, adiabatic surfaces, and hybrid masses. 

With regards to the general features of the adiabatic surfaces, it is possible for spherical bag 
models to reproduce the lattice calculations at small quark separation, but they fail at large R. 
Furthermore, flux tube model or Nambu-Goto string models reproduce the lattice reasonably well 
for intermediate to large quark separations, but do not perform well at small R or in detail at large 
R. 

The flux tube model fits the first excited state, H-^, quite well over a wide range of the quark 
separation (one must ignore the Coulomb term present in Eqn. 25 of Ref. [^). It is, however, 
the only surface to do so at small distance. Furthermore, this may be a fluke due to the particular 
choice of the short distance cutoff of the 7r/i? term employed in Ref. (6^ , 

Of course, n/R splittings are expected to be universal behavior of string systems. Juge, Kuti, 
and Morningstar have carried out a detailed analysis of the relationship of the hybrid surfaces of 
Fig. a to string excitations. They have found that /R splittings are manifest only for very large 
source separation (roughly 4 fermi or greater). This is something of a surprise since one expects a 
phonon-like excitation spectrum on general grounds. 

The lightest flux tube multiplet: (0,1,2) 1 , (0,1,2)^ , and 1++ is larger than that 

observed in lattice computations. Taken together, it appears that the flux tube model is not an 
effective model of gluonic excitations. Alternatively, bag models that have a lowest energy TE mode 
yield multiplets in agreement with the lattice. Yet these models suffer from internal consistency 
issues (center of mass motion, bag distortion, double counting degrees of freedom) that render their 
detailed predictions inaccurate or suspect. 
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Simply achieving the quantum numbers of the adiabatic potentials is not possible in many 
models. Thus, for example, models that employ single spinless gluons do not contain sufficient 
degrees of freedom to reproduce the adiabatic potential spectrum. In particular, transverse “single 
bead” flux tube models (Ref. [^) cannot make lig or states while three dimensional “bead” 
models (Ref. (83| ) cannot make states. Thus including gluon spin is a minimal necessity in 
this class of models (although the level ordering problem must be overcome as discussed above). 
Furthermore, including this spin is not sufficient to obtain the correct level ordering [^, as it 
appears that three-body interactions are required, at least in some models 88 . Alternatively, 
models with constituent gluons with = 1“ in relative S-waves with qq systems fail to describe 
the observed lattice states. However, if the dynamics favors a P-wave coupling then the lowest 
hybrid multiplet can be explained if the qq system is in an S-wave and the heavier lattice multiplet 
can be explained if the quark pair is in a P-wave. If this promise is fulfilled, constituent glue models 
must still strive to incorporate the effects of multiple gluons, which are presumably required to 
obtain string-like large distance behavior [IM] . 


4 Hybrid Decay Models 

Hybrid decay mechanisms were considered from the beginnings of the development of hybrid models. 
Tanimoto categorized these in terms of three possibilities [102| : (i) decays of the type hybrid ^ 
hybrid + meson; these can occur via a mechanism for quark pair production, such as the decay 
model im, [IE; (ii) hybrid ^ hybrid -h meson ^ meson + meson, a two-step process in which 
the virtual hybrid undergoes a gluonic de-excitation to a (conventional) meson state; (hi) hybrid 
^ meson -h meson, where the valence gluon produces a quark pair. 

An important early observation was that a 1 ^ hybrid cannot decay to identical S-wave mesons 
[8^[8^|102| . This result was obtained in the non-relativistic limit employing a decay vertex of the 
form Early computations of this sort found that TM hybrids tended to have typical hadronic 

widths, while TE hybrids tend to be narrow [102 . 


4.1 Flux Tube Decay Model 


The difficulties bag models had in predicting detailed properties of the conventional mesons led to a 
waning of interest in them and an increasing reliance on flux tube models. In fact, shortly after its 
introduction, the flux tube model of meson structure was extended by Isgur, Kokoski, and Patou 
(IKP) to provide a description of meson (m and hybrid [106| decays. The transition operator was 
envisioned as arising due to the quark hopping term of the lattice QCD Hamiltonian. The lowest 
terms in the expansion of this operator are 


Hhop E V’n® • + a E] 


(41) 


n^ji 


n^ji 


where /i is a unit vector along the string. If one assumes a smooth string then the first term 
dominates as the lattice spacing gets small and one has a ^Si strong decay operator. Alternatively, 
if the string is rough then the first term averages to zero upon summing over all local string 
orientations and the second term dominates, yielding a strong decay operator. The authors of 
Ref. [105 , 106| assumed the second scenario since it has a long history of phenomenological success. 

Flux tube degrees of freedom were incorporated by assuming factorization: 
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({.. .}bd- {.. .}bd\0\{.. « (M; MfPolbU^) ■ {{.. .}; {.. .}|{...}). (42) 

Here refers to the collection of (phonon) quantum numbers required to describe the flux 

tube and (9 is a transition operator. The first matrix element on the right hand side is a typical 
^Po mesonic decay overlap. The second represents the overlap of the gluonic/flux tube degrees 
of freedom. Assuming that the quark pair creation occurs at a transverse distance y_i from the 
interquark axis of the parent meson yields the results 


for meson decay and 


({0. 



(43) 

({0.. 



(44) 


for hybrid decay. The factor / is a computable constant of order unity. The extra factor of y± in 
the hybrid decay vertex forces the decay to pairs of identical *S-wave mesons to be zero. This gives 
rise to a selection rule that states that hybrids cannot decay to meson with identical S-wave spatial 
wavefunctions. This is because the string angular momentum cannot be absorbed by the relative 
coordinate between the final state mesons and therefore must go into a qq angular momentum. We 
stress that this result hinges on the assumed factorization of the quark pair production operator 
from the gluonic degrees of freedom. A more general discussion of selection rules in hybrid decays 
is contained in Ref. (I07l 


The observation that flux tube motion can affect the center of mass of the quarks that was 
mentioned in section 3.3.2 was expanded upon by Close and Dudek who commented that radiative 
decays of hybrid mesons can proceed because the recoil of the radiating quark affects the string 
degrees of freedom giving a nonzero overlap of the flux tube wavefunction with the ground state 
flux tube wavefunctions of ordinary mesons |108| . 

A similar scheme involving the emission of pointlike pions may be used to compute hybrid 
decays to final states such as up 109| . The most striking result here is that this decay mechanism 
evades the selection rule discussed above. 


4.2 The PSS Hybrid Decay Model 

An alternative hybrid decay model has been developed in which the gauge field in the interaction 

This model is quite different 


no 
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—g f • Ap; was mapped to phonon degrees of freedom 
from that of IKP because it explicitly correlates quark pair production with gluonic degrees of 
freedom. The resulting decay vertex is given by 


Hint = ^ cos(7re?;“ bl{^rQQ)(T ■ exirgg) (a^x “ (45) 

where the e(r) are polarization vectors orthogonal to r. The qq creation occurs on a line joining the 
original QQ quarks, smeared over the transverse size of the flux tube. The spin operator contracts 
with the flux tube phonon polarization vector, which is absent in the IKP model. Finally, the 
decay amplitude vanishes when the final mesons are identical due to the nodal structure in the 
vector potential. This is true for any single-phonon hybrid in an odd mode. Thus one obtains the 
selection rule: low-lying hybrids do not decay to identical mesons. This subsumes the selection rule 
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of IKP so that none of their qualitative conclusions are changed. However the model predicts, for 
example, that hybrids do not decay to pairs of identical P-wave mesons. 

Another rule, the “spin selection” rule, exists: if the qq in either hybrid or conventional mesons 
are in a net spin singlet configuration then decay into final states consisting only of spin singlet 
states is forbidden. This rule follows because pair creation is spin-triplet and appears to be a 
universal feature in all non-relativistic decay models. 

Notice that the spin selection rule can be used to distinguish conventional and hybrid vector 
mesons. This is because conventional vector mesons have = ^Si or ^Di while hybrid 

vectors have their quarks coupled to a spin singlet. Thus, for example, the decay of hybrid is 
forbidden to tt/ii whereas ttui is allowed; this is the reverse of the case of ^Si conventional mesons 
where the vrai channel is relatively suppressed and irhi is allowed. 

Finally, the IKP and PSS decay models can be compared to the UKQCD lattice results for tti 
decay with the modification introduced by Burns and Close (discussed in section 2.6.3). The model 
widths reported in Table 16 are obtained from Ref. (ng,|m]^ 


mode 

Tpss (MeV) 

Tikp (MeV) 

Tlattice (MeV) 

TTi fin 

10-18 

14 

25 

TTi ^ biTT 

40-78 

51 

80 


Table 16: Comparison of lattice and model hybrid decay modes. 


4.3 Hybrid Photocoupling and Production 

Model hybrid decay computations focus on OZI allowed decays, and thus charmonium or bottomo- 
nium decays are to open flavor final state mesons. In the case of light flavors, an OZI allowed decay 
to a final state such as up will permit an estimate of the radiative transition to Try via the vector 
meson dominance model. The phenomenology of this case was examined in Ref. |112| ; the authors 
noted that hybrid photoproduction could be significant in the case of cj or tt exchange because the 
incoming p is replaced by a photon and because the exchanged particles are off-shell. 

The cross section for deep exclusive electroproduction of an exotic 1 ^ hybrid was estimated 
in the Bjorken regime by Anikin et al. [HI- This cross section was found to be large, and in fact 
scales like as is usual for meson electroproduction. The authors also noted that forward- 

backward asymmetry in the production of tt and p mesons can serve as a useful signal for hybrid 
meson production. 

We are aware of one attempt to compute the hidden flavor decay of hybrid mesons (m. The 
model assumed hybrids constructed with constituent gluons such that the TE mode is the lowest 
lying. The electromagnetic transition was modelled in bound state perturbation theory with the 
valence gluon being absorbed by quarks or by the Coulomb interaction with photon emission from 
a quark line. The authors were not able to predict absolute rates so that comparison to the lattice 
results of sect ion [2. 6. 4| is not possible. The rates for all the low lying TE hybrid multiplet transitions 
were found to be comparable. 

The utility of polarized photoproduction of hybrid mesons was analyzed in Ref. [HI- The 
authors examined the process yp ^ p^n^n at low momentum transfer via pion exchange. They 
assumed a simple form for the photon-pion-hybrid vertex and fixed the ratio of partial waves 
according to the vector meson dominance model. The computed cross sections indicated that 
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corrections due to absorption are required and that linear photon polarization permits separating 
the exotic wave, even without partial wave analysis. 


5 Experimental Situation 


Experimental data on exotic hybrid mesons goes back nearly three decades, and continues until 
today. An extensive review on spectroscopy in general can be found in reference |116| , while a 
review focussing only on hybrid mesons is given in reference The experimental situation 

has evolved since these previous reviews, with new data from several experiments and detailed 
systematic studies of production mechanisms made possible by very-large data sets. To date, there 
is experimental evidence for three exotic quantum number states [118| . All have the same quantum 
numbers, the 7ri(1400), 7ri(1600) and 7ri(2015), and there are experimental and 

interpret at ional issues surrounding all of them. The field is now on the cusp of a new experiment 
to study the photoproduction of hybrids, and a few years away from a new antiproton experiment. 
Thus it seems timely to take stock of the current situation. 

In order to search for exotic hybrids, it is first necessary to catalog the possible decay modes. 


These are given in Table 17 where the widths come from model calculations (see Section 4.2). 
The decay modes are obtained by simple quantum number counting arguments. The list is not 
exhaustive, but covers the most likely and most detectable modes. 


Name 


Total Width MeV 

Allowed Decay Modes 



PSS 

IKP 


TTl 

1-+ 

81 - 168 

117 

biTT, np, tt/i, Trr/, irr]', pai, 7rr?(1295) 

m 

1-+ 

59 - 158 

107 

vrai, 7ra2, ??/i, ??/ 2 , 7r7r(1300), rirj', KKf, KKf 


1-+ 

95 - 216 

172 

KK^, KK^, KK*, 7]7]' 

bo 

0+- 

247 - 429 

665 

7r7r(1300), irhi, pfi, pbi 

ho 

0+- 

59 - 262 

94 

irbi, rjhi, it'it'(1460) 

h'o 

0+- 

259 - 490 

426 

KK{U60), KK^, rjhi 

b2 

2+- 

5-11 

248 

nai, 7ra2, nhi, pp, pbi, pfi 

h2 

2+- 

4-12 

166 

7rp, TT^I, TjUJ^ (jjbi 

h'2 

2+- 

5-18 

79 

KK^, KK^, KK^, phi 


Table 17: Exotic quantum number hybrid widths and possible decay modes 


5.1 Exotic Quantum Number states 


Data on exotic-quantum-number mesons have come from both diffractive production using incident 
pion beams, antiproton annihilation on protons and neutrons and charmonium decays. Diffractive 
production is schematically shown in Fig. 11 A pion beam is incident on a proton (or nuclear) 
target, which recoils following a t-channel exchange. The process can be written down in the 
reflectivity basis m where the production factorizes into two non-interfering amplitudes—positive 
reflectivity (e = -h) and negative reflectivity {e — The absolute value of the spin projection 
along the z-axis is M, and is taken to be either 0 or 1 (it is usually assumed that contributions from 
M larger than 1 are small and can be ignored |120| ). It can be shown in this process that naturality 
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of the exchanged particle can be determined by e. Natural parity exchange (n.p.e.) corresponds 
to J^s of 0+, 1“, 2+, • • •, while unnatural parity exchange (n.p.e.) corresponds to of 0“, 1+, 
2“, • • •. For a state which is observed in more than one decay mode, one would expect that the 

TT (beam) X(J^M^) 

\ s = + : n.p.e. 

\ £ = - : u.p.e. 



p (target) P (recoil) 

Figure 11: The diffractive process with an incident pion (tt beam) and a proton (p) target. The 
exchange has z-component on angular momentum M and reflectivity e. The recoil proton and a 
state X of given constitute the final state. For positive reflectivity, the t-channel is a natural 
parity exchange, while for negative reflectivity, it is unnatural parity exchange. 

production mechanism (M^) would be the same for all decay modes. If not, this could be indicative 
of more than one state being observed, or possible analysis problems. 

In antiproton-nucleon annihilation, there are a number of differences between various annihila¬ 
tion processes. For the case of pp, the initial state is a mixture of isospin / = 0 and / = 1. For pn 
annihilation, the initial state is pure I — 1. For annihilation at rest on protons, the initial state is 
dominated by atomic S-waves. In particular, ^Sq and ^Si atomic states, which have = 0 ^ 
and 1 respectively (with a small admixture of P states). For annihilation in flight, the number 
of initial states is much larger and it may no longer make sense to try to parametrize the initial 
system in terms of atomic pp states. 

The combination of initial isospin and final state particles may lead to additional selection rules 
that restrict the allowed initial states. In the case of pp pyr^Tr^, the annihilation is dominated 
by ^Sq initial states (J^^ = 0 ^). For the case of pn pyr^Tr”, quantum numbers restrict this 
annihilation to occur from the ^Si initial states (J^^ = 1 ). In addition, the neutron is bound in 

deuterium, where the Fermi motion introduces substantial p-wave annihilation. Thus, one may see 
quite different final states from the two apparently similar reactions. The results to date came from 
the Crystal Barrel experiment at CERN. New results are expected from the PANDA experiment 
at FAIR sometime after 2020. 

A limited number of results on the photoproduction of exotic mesons have been published, with 
no reported signals to date. Extensive new results are expected with the advent of the GlueX 
experiment at Jefferson Lab in 2015, and the CLAS-12 experiment two years later. 

The decays of Xc niesons should also be a good place to search for hybrid mesons. Results 
from CLEO-c report interesting results in the pV system. More extensive results with much higher 
statistics should be available from BES-III. 

5.2 The yri(1400) 

5.2.1 GAMS results on the yri(1400) 

The GAMS experiment made the first reported observation of an exotic quantum number state, 
using a 40 GeV/c yr“ beam to study the reaction 7T~p pp7T~. They reported a — 1 ^ state 
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Figure 12: The results of a partial-wave analysis of the rj'K~ final state from VES. (a) shows the 
intensity in the partial wave, (b) shows the intensity in the 1 ^ partial wave and (c) shows 
the relative phase between the waves. (Figure reproduced from reference 126| .) 


(the M(1405)) in the r]7T system 121 . They reported a mass of 1.405 ± 0.020 GeV and a width 
of 0.18 ± 0.02 GeV for this state. In the neutral channel, an earlier search found no evidence 


for an exotic quantum-number state [122 


5.2.2 KEK results on the 7ri(1400) 

Somewhat later, an experiment at KEK using a 6.3 GeV/c tt 
r]7T~ system with a mass of 1.3431 ± 0.0046 GeV and a width of 0.1432 ± 0.0125 GeV [123 


beam observed a 1 


state in the 
Since 


this mass and width is very close to that of the strong signal for the a2(1320), the possibility of 
signal leakage into the weak exotic wave from the strong a 2 signal often comes up. 


5.2.3 VES results on the 7ri(1400) 

The VES experiment studied interactions using a 37 GeV/c pion beam. They reported intensity in 
the 1 ^ ri7T~ exotic wave and saw rapid phase motion between the 1 ^ wave and the a2(1320) |124 


12 


as shown in Fig. 
not in the 0“ and 1 


The exotic wave was present in the — 1+ (natural parity) exchange, but 
(unnatural parity) exchange. The observed data could be fit using — 1 ^ 


intensity and the phase motion with respect to the a2(1320) using a Breit-Wigner distribution with 
a mass of 1.316 ± 0.012 GeV and width of 0.287 ± 0.025 GeV. However, VES stopped short of 
claiming an exotic resonance, as they could not unambiguously establish the nature of the exotic 
In a later analysis of the t^tt^ system, they claimed that the peak near 1.4 GeV can be 


wave 


125 


understood without requiring an exotic quantum number meson [126] . 


5.2.4 E852 results on the 7ri(1400) 


prjiT . They 
The state was only produced 


The E852 collaboration used an 18 GeV/c beam to study the reaction 7T~p 
reported the observation of a 1 ^ state in the p'K~ system 12 
in natural parity exchange (M^ = 1^). Fitting to a Breit-Wigner distribution, they found a 
mass of 1.37 ± 0.016^0 030 GeV and a width of 0.385 ± 0.040^0 105 GeV. The intensity and phase 
difference plots are shown in Fig. While the observed exotic signal was only a few percent of the 


dominant a2(1320) strength, they noted that its interference with the a 2 provided clear evidence 
of the existence of the 7ri(1400). This is most clearly shown in plot (d) of the figure where curve 1 
is the phase of the a 2 , curve 2 is the phase of the tti, 3 is an assumed flat background phase, and 
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4 is the resulting phase difference. When their intensity and phase-difference plots were compared 
with those from VES m, they were found to be identical. 



Figure 13: Data from the E852 experiment showing the 7ri(1400) |127| . In (a) is shown the intensity 
of the = 2++ partial wave, dominated by the a2(1320), as a function of rjir mass. In (b) is 
shown the intensity of the exotic 1 ^ wave as a function of the titt mass. Frame (c) shows the 
phase difference between the two partial waves, 2++ and 1 Frame (d) shows the evidence for 
the resonant nature of the 7ri(1400) by looking at the three components that build the total phase 
difference. (Figure reproduced with permission from reference (m|.) 


Due to disagreements over the interpretation of the 1 ^ signal, the E852 collaboration split into 
two groups. The majority of the collaboration published the resonance interpretation, tti (1400) [127 


while a subset of the collaboration did not sign the paper. As this latter group, centered at Indiana 
University, continued to analyze data collected by E852, we will refer to their publications as E852- 
lU to distinguish the work of the two groups. An analysis by the E852-IU group of data for the 
reaction 7T~p npir^ found evidence for the exotic 1 ^ partial wave, but was unable to describe it 


as a Breit-Wigner-like 7ri(1400) resonance 128 . However, a later analysis by the E852 collab¬ 


oration of the same final state and data confirmed their earlier observation of the 7ri(1400) [129 


E852 found a mass of 1.257 ± 0.020 ± 0.025 GeV and a width of 0.354 ± 0.064 ± 0.058 GeV with 
the 7ri(1400) produced via natural parity exchange (M^ = 1+). Much of the discrepancy be¬ 
tween these two works arose from the treatment of backgrounds. The E852 collaboration consider 
no background phase, and attribute all phase motion to resonances. The E852-IU group allow for 
non-resonant interactions in the exotic channel, these background processes are sufficient to explain 
the observed phase motion. 


5.2.5 Crystal Barrel results on the 7ri(1400) 

The Crystal Barrel Experiment studied antiproton-neutron annihilation at rest in the reaction 


pn 


rjTT TT^ 


130 . The Dalitz plot for this final state is shown in Fig. 14 where bands for the 


a2(1320) and p(770) are clearly seen. They reported an exotic = 1 + state, the 7ri(1400), with 
a mass of 1.40 ± 0.020 ± 0.020 GeV and a width of 0.310 ± 0.050^0 030 GeV. While the signal is not 
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obvious in the Dalitz plot, if one compares the difference between a fit to the data without and with 

. While 
Crystal 


the 7ri(1400), a clear discrepancy is seen when the 7ri(1400) is not included (see Fig. 15). While 
the 7ri(1400) was only a small fraction of the a2(1320) in the E852 measurement 1 
Barrel observed the two states produced with comparable strength. The Crystal Barrel group 
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Figure 14: (Color on line.) 


rfTT TT^ from the Crystal Barrel experiment 130 


The Dalitz plot of m?{r]7r^) versus ) for the reaction pn 

The bands for the a2(1320) are clearly seen in 


both r]7v^ and piv , while the p(770) is seen in the tt^tt invariant mass. 


also looked at the reaction pp 


rjTV^TV^ 


131 


Here, a weak signal was observed for the 7ri(1400) 


(relative to the a2(1320)) with a mass of 1.360 ± 0.025 GeV and a width of 0.220 ± 0.090 GeV. In 
I = 0 pp annihilations, the a2(1320) is produced strongly from the ^Sq atomic state. However, the 
pn is isospin 1 and state is forbidden. Thus, the strong a 2 production from pp is suppressed in 
pd annihilations—making the 7ri(1400) production appear enhanced relative to the a2(1320) in the 
latter reaction. 


5.2.6 OBELIX results on the 7ri(1400) 


Both the OBELIX [132 experiment reported on the observation of the 7ri(1400) in the reaction 
pp Att. They both claimed the observation of the 7ri(1400) decaying to piv final states, however 
there is some concern about the production mechanism. This same result was reported by the 
Crystal Barrel Experiment [133] in conference only. The pir signal arises from annihilation from 
the p-wave initial state, while the signal in pn come from the initial state. Thus, it is unlikely 
that the exotic state seen in rjiT and that seen in pn are the same. The origin of these may not be 
due to an exotic resonance, but rather some re-scattering mechanism that has not been properly 
accounted for. 


5.2.7 COMPASS results on the 7ri(1400) 

The COMPASS experiment used a 191 GeV/c n~ beam to study the production of both the pn 
and p'n final states 134 . They note that in the pn system, the exotic 1 ^ wave “shows a compact 
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Figure 15: (Color on line.) The difference between the fit and the data in the D alitz plot of m^{r]7T^ 
versus t/tt” for the reaction pn r]7T~7T^ from the Crystal Barrel experiment 130 . (a) Does not 


include the 7ri(1400) while (b) does include the 7ri(1400). There are clear systematic discrepancies 
present in (a) that are not present when the 7ri(1400) is included. 


peak of 400 MeV/c^ width, centered at a mass of 1.4 GeV/c^”. 
resonant nature of this structure. 


They make no comments on the 


5.2.8 Interpretation of the 7ri(1400) 

Interpretation of the 7ri(1400) has been problematic. Its mass is lower than most predicted values for 
hybrid mesons, and its observation in only a single decay mode (ijn) is not consistent with models of 
hybrid decays. Donnachie and Page showed that the 7ri(1400) could be an artifact of the production 
dynamics. They demonstrated that is possible to understand the 7ri(1400) peak as a consequence 


of the 7ri(1600) (see Section 5.3) interfering with a non-resonant Deck-type background with an 
appropriate relative phase 135 . Zhang [136| considered a molecular picture where the 7ri(1400) 
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was an 77(1295)7r molecule. However, the predicted decays were inconsistent with the observations 
of the 7ri(1400). 

Szczepaniak et al Eg considered a model in which t-channel forces could give rise to a back¬ 
ground amplitude which could be responsible for the observed 7ri(1400). In their model, meson- 
meson interactions which respected chiral symmetry were used to construct the t/tt p-wave interac¬ 
tion much like the tttt 5 -wave interaction gives rise to the a meson. They claimed that the 7ri(1400) 
was not a QCD bound state, but rather dynamically generated by meson exchange forces. 

Close and Lipkin noted that because the SU(3) multiplets to which a hybrid and a multiquark 
state belong are different, rjir and decays might be a good way to distinguish them. They found 
that for a multiquark state, the titt decay should be larger than 77 ^, while the reverse is true for a 
hybrid meson |138| . A similar observation was made by Chung 139| who noted that in the limit of 
the 77 being a pure octet state, the decay of an octet 1 ^ state to an titt p-wave is forbidden. Such 
a decay can only come from a decuplet state. Given that the pseudoscalar mixing angle for the p 
and p' are close to this assumption, they argue that the 7ri(1400) is qqqq in nature. 

While the interpretation of the 7ri(1400) is not clear, most analyses agree that there is intensity 
in the 1 ^ wave near this mass. A summary of all reported masses and widths for the 7ri(1400) 
is given in Table All are reasonably consistent, and even the null observations of VES and 


E852-IU all concur that there is strength near 1.4 GeV in the exotic wave. However, the E852 
and VES results can be explained as either non-resonant background pTl , or non-resonant Deck 
amplitudes |135| . Another possibility is the opening of meson-meson thresholds, such as /i(1285)7r. 
Unfortunately, no comparisons of these hypothesis have been made with the pN data (owing to the 
lack of general availability of the data sets), so it is not possible to conclude that they would also 
explain those data. 


Mode 

Mass (GeV) 

Width (GeV) 

Experiment Re 

hrence 

p7r“ 

1.405 ±0.020 

0.18 ± 0.02 

GAMS 

121 


p7r“ 

1.343 ±0.0046 

0.1432 ±0.0125 

KEK 

123 


p7r“ 

1.37 ±0.016 

0.385 ± 0.040 

E852 

127 


pTT^ 

1.257 ±0.020 

0.354 ± 0.064 

E852 

129 


pTT 

1.40 ± 0.020 

0.310 ±0.050 

CBAR 

130 


pvr^ 

1.36 ±0.025 

0.220 ±0.090 

CBAR 

131 


piT 

1.384 ±0.028 

0.378 ± 0.058 

Obelix 

132 


pTT 

~ 1.4 

~ 0.4 

CBAR 

133 


pTT 

1.354 ± 0.025 

0.330 ± 0.035 

PDG 

118 



Table 18: Reported masses and widths of the 7ri(1400) from the GAMS, KEK, E852, Crystal 
Barrel (CBAR) and Obelix experiment. Also reported is the 2014 PDG average for the state. Not 
reported here is the COMPASS observation of intensity in the 1 ^ partial wave near a mass of 
1400 MeV p4| . 

5.3 The 7ri(1600) 

While the issue of only one decay mode (p 7 r) and the low mass have made interpretation of 7ri(1400) 
nature difficult, a second = 1 ^ state has a much richer set of observations, both in production 

and decay, and in statistics. This other state, the 7ri(1600), has been observed in diffractive 
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production using incident 7 r“ beams where its mass and width have been reasonably stable over 
several experiments and observed decay modes. There is some evidence in pp annihilation as well 
as Xci decays. These positive results have been reported from VES, E852, COMPASS, Crystal 
Barrel, CLEO-c and others, and are discussed below. Negative results have been reported in 
photoproduction experiments by CL AS. 


5.3.1 VES Results on the 7ri(1600) 

In addition to their study of the r]7v~ system, the VES collaboration also examined the ryV” system. 
Here they observed a = 1 ^ partial wave with intensity peaking at a higher mass than the 
7ri(1400) [m . However, as with the rjiT system, they did not claim the discovery of an exotic- 
quantum-number resonance. VES later reported a combined study of the 77 ^“, fi7r~ and p^7v~ 
final states [l40| , and reported a “resonance-like structure” with a mass of 1.62 ± 0.02 GeV and 
a width of 0.24 ± 0.05 GeV decaying into p^7T~. They also noted that the wave with — 1 ^ 
dominates in the 77 ^“ final state, peaking near 1.6 GeV and observed a small 1 ^ signal in the 
fi7v~ final state. 

VES also reported on the uj7t~7v^ final state fliTj |l42| . In a combined analysis of the rj^n , 
biTT and p^7v~ final states, they reported the 7ri(1600) state with a mass of 1.61 ± 0.02 GeV and a 
width of 0.29 ± 0.03 GeV that was consistent with all three final states. To the extent that they 
observed these states, they also observed all three final state produced in natural parity exchange 
(M^ = 1+). They were also able to report relative branching ratios for the three final states as 
given in equation 1 ^ 


biTT : piT \ pTV \ = 1 : 1 ± 0.3 : 1.5 ± 0.5. 


(46) 


However, there were some issues with the pir final state. Rather than limiting the rank of the 
density matrix as was done in [143 , 144] , they did not limit it. This allowed for a more general fit 
that might be less sensitive to acceptance affects. In this model, they did not observe any significant 
structure in the 1 ^ pn partial wave above 1.4 GeV. However, by looking at how other resonances 
were produced, they were able to isolate a coherent part of the density matrix from which they 
found a statistically significant 1 ^ partial wave peaking near 1.6 GeV. While VES was extremely 
careful not to claim the existence of the pn decay of the the 7ri(1600), in the case that it exists. 


they were able to obtain the rates given in equation 46 



Figure 16: The results of a partial wave analysis on the 77 ^“ final state from VES. (a) shows 
he 2 ++ partial wave in cjp, (b) shows the 1 ^ partial wave in biir and (c) shows the interference 

between them. (Figure reproduced with permission from reference (H.) 

The VES results have been summarized in a review of all their work on hybrid mesons [126 . 
This included an updated summary of the 7ri(1600) in all four final states, p'lr p 7 r, biir and /itt. In 
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Figure 17: The results of a partial wave analysis on the bin final state from VES. (a) shows he 
partial wave, (b) shows the 1 ^ partial wave and (c) shows the interference between them. (Figure 
reproduced with permission from reference 126| .) 
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Figure 18: The results of a partial wave analysis on the fin final state from VES. (a) shows he 1++ 
partial wave, (b) shows the 1 ^ partial wave and (c) shows the interference between them. (Figure 
reproduced with permission from reference 126| .) 


the r]'n final state (Fig. 16), they note that the 1 ^ wave is dominant. While they were concerned 
about the nature of the higher-mass part of the 2++ spectrum (a2(1700) or background) they find 
that a resonant description of 7ri(1600) was possible in both cases. For the case of the bin final 
state (Fig. 0 , they found that the contribution of a 7ri(1600) resonance is required. In a combined 
fit to both the rj'n and bin data, they find a mass of 1.56±0.06 GeV and a width of 0.34±0.06 GeV 
for the 7ri(1600). In the fin final state (Fig. [^, they find a resonant description of the 7ri(1600) 
with a mass of 1.64 ± 0.03 GeV and a width of 0.24 ± 0.06 GeV which they note is compatible 
with their measurement in the previous two final states. They also note, that in contradiction with 
E852 146 , they find no significant 1 ^ intensity above a mass of 1.9 GeV (see Section 5.4). For 


the pn final state, they are unable to conclude that the 7ri(1600) is present. 

They note that the partial-wave analysis of the n^n~n~ system finds a significant contribution 
from the = 1^ wave in the pn channel (2 to 3% of the total intensity). Some of the models 
in the partial-wave analysis of the exotic wave lead to the appearance of a peak near a mass of 
1.6 GeV which resembles the 7ri(1600). However, the dependence of the size of this peak on the 
model used is significant 142 . They note that because the significance of the wave depends very 


strongly on the assumptions of coherence used in the analysis, the results for 37r final states on the 
nature of the 7ri(1600) are not reliable. 

To obtain a limit on the branching fraction of 7ri(1600) decay to p7r, they looked at their results 
of the production of the 7ri(1600) in the charge-exchange reaction to rj'n^ versus that of the ri'n~ 
final state. They conclude that the presence of the 7ri(1600) in ri'n~ and its absence in rj'n^ preclude 
the formation of the 7ri(1600) by p exchange. From this, they obtain the relative branching ratios 


for the 7ri(1600) as given in equation 47 


bin : fin : pn : pn = 1.0 ± .3 : 1.1 ± .3 :< .3 : 1. 


(47) 


40 


































5.3.2 E852 Results on the 7ri(1600) 


Using an 18 GeV/c tt 
partial wave analysis of the final state 

intermediate states and observed a = U 


143 


beam incident on a proton target, the E852 collaboration carried out a 

They saw both the p^7r~ and /2(1270)7r“ 
state which decayed to p7r, the 7ri(1600). The 
7ri(1600) was produced in both natural and unnatural parity exchange ( — 1+ and — 0“, 1“) 


with apparent similar strengths in all three exchange mechanisms (see Fig. 19). In Ref. (m. they 
noted that there were issues with the unnatural exchange production. The signal in the = 1“ 
wave exhibited very strong model dependence and nearly vanished when larger numbers of partial 
waves were included. The signal in the — 0“ partial wave was stable, but its peak was above 
1.7 GeV. They noted that the unnatural-parity exchange is expected to die off at higher energies, 
so their results are not at odds with those of VES, where natural parity exchange dominates. In 
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Figure 19: The production of the 1 partial wave as seen in the 7r'’"7r“7r“ final state by E852. (a) 
shows the unnatural parity exchange (M^ = 0“,1“) while (b) shows the natural parity exchange 
(M^ = 1^). (Figure reproduced with permission from reference (i^.) 

the E852 data, the unnatural parity exchange waves make up a small fraction of the total signal. In 
unnatural parity exchange, they found no significant waves, which made a study of phase motion of 
the 1 ^ in this sector problematic. Thus, in their analysis, they only considered the natural parity 
exchange. There, they found the 7ri(1600) to have a mass of 1.593±0.08^ 

0.168 ± 0.020tg;J^^ GeV. In Fig, 


20 


are shown the intensity of the 1“ 


^0 047 GeV and a width of 
and 2 ^ (772(1670)) partial 


waves as well as their phase difference. The phase difference can be reproduced by two interfering 
Breit-Wigner distributions and a flat background. 

In a follow-up analysis, E852 also studied the reaction 7r“p ^ pryV” to examine the ri^7T~ 
final state They observed, consistent with VES (T^ , that the dominant signal was the 

1 ^ exotic wave produced dominantly in the — 1^ channel, implying only natural parity 

exchange. They found the signal to be consistent with a resonance, the 7ri(1600) and found a mass 
of 1.597 iO.OlOlg^ GeV and a width of 0.340 ±0.040 ±0.050 GeV. The results of the E852 PWA 
are shown in Fig. 


21 


where the wave is the 1 the Dj^ corresponds to the 2++ a 2 and the G+ 
corresponds to the a 4 . Clear phase motion is observed between both the 2++ and 4++ wave 
and the 1 ^ and the 4++ wave. 

E852 also looked for the decays of the 7ri(1600) to 5i7r and /itt. The latter was studied in the 
reaction pr]7T^7T~7T~ with the fi being reconstructed in its 7777 + 77 “ decay mode The 

771 (1600) was seen via interference with both the 1++ and 2 ^ partial waves. It was produced via 
natural parity exchange (M^ = 1+) and found to have a mass of 1.709 ± 0.024 ± 0.041 GeV and 
a width of 0.403 ± 0.080 ±0.115 GeV. A second 771 state was also observed in this reaction (see 


Section 5.4). 
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Figure 20: The results of a PWA to the 7t^7t~7t~ final state from E852. (a) shows the intensity 
of the = 1 ^ wave, (b) shows the 2 ^ and (c) shows the phase difference between the two. 

The solid curves are fits to two interfering Breit-Wigner distributions. In (d) are shown the phases 
of the two Breit-Wigner distributions and (1,2) and a flat background phase (3) that combine to 
make the curve in (c). (Figure reproduced with permission from reference mi.) 



IVl(Ti'jt-) (GeV) 

Figure 21: Results from E852 on the 77 ^“ final state, (a) shows the 1 ^ partial wave, (b) shows 
the partial wave (an and (c) shows the phase difference between these, (d) shows the 2 ++ 
partial wave (a2(1320)), while (e) shows the a 4 and (f) is the phase difference. (Figure reproduced 
with permission from reference |145| .) 
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The biTT final state was studied by looking at the reaction tt p ^ ujtt ttV? with the bi re- 

The 7ri(1600) was seen interfering with the and 


constructed in its ujtt decay mode 147 


partial waves. In 6i7r, they measured a mass of 1.664 ± 0.008 ± 0.010 GeV and a width of 
0.185 ± 0.025 ± 0.028 GeV for the 7ri(1600). However, the production mechanism was seen to be 
a mixture of both natural and unnatural parity exchange, with the unnatural being stronger. As 
with the /iTT, they also observed a second tti state decaying to biTv (see Section 5.4). 


final state 

production (M^) 

dominant 

piT 

0-, 1-,1+ 

npe ^ upe 

p'tt 

1+ 

npe 

flTT 

1+ 

npe 

b\TT 

0-,1-,1+ 

upe > npe 


Table 19: The production mechanisms for the 7ri(1600) as seen in the E852 experiment. Also shown 
is whether the natural parity exchange (npe) or the unnatural parity exchange (upe) is stronger. 


The fact that E852 observed the 7ri(1600) produced in different production mechanisms, de¬ 
pending on the final state, is somewhat confusing. A summary of the observed mechanisms is given 
in Table In order to understand the variations in production, there either needs to be two 
nearly-degenerate 7ri(1600)s, or there is some unaccounted-for systematic problem in some of the 
analyses. 



Figure 22: (Color on line) The E852-IU PWA solutions for the 1 ^ partial wave for the 

final state (a) and its interference with the 2++ partial wave (b). See text for an explanation of the 

labels. (Figure reproduced with permission from reference [148| .) 


The E852-IU group analyzed an E852 data set that was an order of magnitude larger than 
that used by E852 in Refs. 


143 


144 


__ In this larger data set, they looked at the reactions 7T~p 

n7T^7T~7T~ and 7T~p nTT^TT^TT^ and carried out a partial wave analysis for both the tt^tt~tt~ and 
the TT~TT^TT^ final states. This yielded solutions that were consistent with both final states 148 . In 


this analysis, they carried out a systematic study of which partial waves were important in the fit. 
When they used the same wave set as in the E852 analysis [143 , 144 , they found the same solution 
showing a signal for the 7ri(1600) in both final states. However, when they allowed for more partial 
waves, they found that the signal for the 7ri(1600) vanished. Fig. 22 shows these results for the 
TT~TT^TT^ final state, while Fig. shows the results for the tt^tt~tt~ final state. In both figures, the 
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Figure 23: (Color on line) The E852-IU PWA solutions for the 1 ^ partial wave for the tt^tt^tt” 
final state (a) and its interference with the 2++ partial wave (b). See text for an explanation of the 
labels. (Figure reproduced with permission from reference |148| .) 


“low wave” solution matches that from E852, while their “high wave” solution shows no intensity 
for the 7ri(1600) in either channel. An important point is that in both their high-wave and low-wave 
analyses, the phase difference between the exotic 1 ^ wave and the 2^^ wave are the same (and 
thus the same as in the E852 analysis). While not shown here, the same is also true for the 1 ^ 
and 2 ^ waves. 


7r2(1670) 


= 0+ 


= 1+ 

M" = 1- 

Decay 

L 

H 

L 

H 

L H 

(/2vr)s 

X 

X 

X 

X 

X 

(/27r)D 

X 

X 

X 

X 


[(TrTr)^]^ 

X 

X 


X 


{PT^)P 

X 

X 


X 


{pp)f 


X 


X 


(Md 


X 


X 



Table 20: The included decays of the 7r2(1670) in two analyses of the Stt final state. “L” is the 
wave set used in the E852 analysis (l43l[^. “H” is the wave set used in the higher statistics 


analysis 148 


E852-IU carried out a study to determine which of the additional waves in their “high wave” set 
were absorbing the intensity of the 7ri(1600). They found that the majority of this was due to the 
inclusion of the pn decay of the 7r2(1670). The partial waves associated with the 7r2(1670) in both 
analyses are listed in Table While the production from — 0^ is similar for both analyses, 
the E852 analysis only included the 7r2(1670) decaying to /27r in the — 1+ production. The 
high-statistics analysis included both /27r and pir in both production mechanisms. The PDG 118| 


lists the two main decays of the 7r2(1670) as /27r (56%) and pir (31%), so it seems odd to not include 
this latter decay in an analysis including the 7r2(1670). Fig. shows the results of removing the 
pn decay from the “high wave” set for the 7t^7t~7t~ final state. This decay absorbs a significant 
portion of the 7ri(1600) partial wave. 

In the E852-IU analyses, the fact that the phase motion of the 1 ^ exotic wave relative to 
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Figure 24: Results from the E852-IU analysis showing the 1 ^ intensity for the charged mode for 
the high-wave set (filled circles), the modified high-wave set (filled squares), and the low-wave set 
(open circles). In the modified high-wave set the two pn decays of the 7r2(1670) were removed from 
the fit. (Figure reproduced with permission from reference [148| .) 


other partial waves agrees with those differences as measured by E852, and are the same in both 
the high-wave and low-wave analyses is intriguing. This could be interpreted as a 7ri(1600) state 
which is simply absorbed by the stronger 7r2(1670) as more partial waves are added. However, 
given the small actual phase difference between the 1 ^ and 2 ^ partial waves (see Fig. 20), the 
opposite conclusion is also possible, particularly if some small non-zero background phase were 
present. Here, the 1 ^ signal is due to leakage from the stronger 7^2 and no 7ri(1600) is needed in 
the piT final state. 

While the results on pn between E852 and VES seem at odds, we believe that these discrepancies 
are the result of the assumptions made in the analyses. These assumptions then manifest themselves 
in the interpretation of the results. The VES analyses fit both the real and imaginary parts of their 
amplitudes independently. However, for analytic functions, the two parts are not independent. 
Not using these constraints can lead to results that may be unphysical, and at a minimum, are 
discarding important constraints on the amplitudes. In E852, many of their results rely on the 
assumption of a flat background phase, but there are many examples where this is not true. Thus, 
their results are biased toward a purely resonant description of the data, rather than a combination 
of resonant and non-resonant parts. It is also somewhat disappointing that E852 is unable to make 
statements about relative decay rates, or carry out a coupled channel analysis of their many data 
sets. Our understanding is that this is due to issues in modelling the rather tight trigger used in 
collecting their data. 


5.3.3 Crystal Barrel results on the 7ri(1600) 


An analysis of Crystal Barrel data at rest for the reaction pp u7t^7t~7t^ was carried by some 
members of the collaboration [149 . They reported evidence for the 7ri(1600) decaying to biir from 
both the ^Sq and ^Si initial states, with the signal being stronger from the former. The total signal 
including both initial states, as well as decays with 0 and 2 units of angular momentum accounted 
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for less than 10% of the total reaction channel. The mass and width were found consistent (within 
large errors) of the PDG value, and only results with the mass and width fixed to the PDG values 
were reported. Accounting for the large rate of annihilation to of 13%, this would imply 

that pp 7ri(1600)7r accounts for several percent of all pp annihilations. 


5.3.4 CLEO-c results on the 7ri(1600) 


In CLEO-c, the decays of the Xci to both p7T^7T~ and t/'tttt were studied. An amplitude analysis of 
the p^TTTT final state showed a 4(j signal for an exotic p-wave in the p'n system. Results from this 
analysis are shown in Fig. While the analysis could not conclude if the P-wave is resonant, if 
it is fit by a Breit-Wigner amplitude, the data can be well described by a 7ri(1600) with a mass of 
1670 ± 30 ± 20 MeV/c^ and a width of 240 ± 50 ± 60 MeV/c^ 


150 



Figure 25: (Color on line.) Results from CLEO-c for the decay of the Xci to The figure 

shows the fit intensity in 77 V system. The solid (red) fit curve shows the intensity in the exotic 
1 ^ partial wave. There is clear peaking near the mass of the 7ri(1600). (Figure reproduced with 
permission from reference [150| .) 


5.3.5 CL AS results on the 7ri(1600) 

The CLAS experiment at Jefferson Lab studied the reaction 779 ^ 7 r+ 7 r+ 7 r“(n )^^55 to look for the 
production of the 7ri(1600) mi. The photons were produced by bremsstrahlung from a 5.7 GeV 
electron beam. While there were significant contributions from baryon resonances in their data, 
they attempted to remove this by selective cuts on various kinematic regions. The results of their 
partial-wave analysis show clear signals for the ai(1270), the a2(1320) and the 7r2(1670), but show 
no evidence for the 7ri(1600) decaying into three pions. They place an upper limit of the production 
and subsequent decay of the 7ri(1600) to be less than 2% of the a2(1320). Their results imply (i) 
the 7ri(1600) is not strongly produced in photoproduction, (ii) the 7ri(1600) does not decay to 37 r, 
or (hi) both. Using a much larger data sample, CLAS has continued to search for the 7ri(1600) in 
photoproduction through the 37 r decay mode, however the preliminary results are still consistent 
with the earlier work |152| . They propose that either there is a preference for production through 
mechanisms other than pion exchange, or that the 7ri(1600) does not decay to 37 r. 
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Figure 26: The results from CLAS of a partial-wave analysis photoproduction data of the Stt final 
state. Intensity is seen in the 2 ^ partial wave, (a) and (b), as well as the partial wave (c). 
In the 1 ^ exotic wave, (d), no intensity is observed. (Figure reproduced with permission from 
reference PD) 


5.3.6 COMPASS results on the 7ri(1600) 


The COMPASS experiment has reported their first study of the diffractively produced Stt final 
state |153 , 154 . They used a 190 GeV/c beam of pions to study the reaction 7 T~Pb 7 t~7t~7t^X. 
In their partial-wave analysis of the Stt final state, they observed the 7ri(1600) with a mass of 
1.660 ± 0 . 010^0 004 GeV and a width of 0.269 ± O.O211 qq 04 GeV. The 7ri(1600) was produced 
dominantly in natural parity exchange (M^ = 1 ^) although unnatural parity exchange also seemed 
to be required. However, the level was not reported. The wave set (in reference 
to be somewhat larger than that used in the high-statistics study of E852-IU 


154 


148 


used appears 
Thus, in the 

COMPASS analysis, the pn decay of the 7r2(1670) does not appear to absorb the exotic intensity in 
their analysis. They also report on varying the rank of the fit with the 7ri(1600), with the results 
being robust against these changes. One point of small concern is that the mass and width that 
they extract for the 7ri(1600) are essentially identical to those for the 7r2(1670). For the latter, they 
observed a mass of 1.658 ± 0 . 002^0 oog GeV and a width of 0.271 ± 0.009^0 024 GeV. However, the 
strength of the exotic wave appears to be about 20 % of the 7 r 2 , thus feed through seems unlikely. 
Results from their partial-wave analysis are shown in Figs. 27 and 28, These show the 1 ^ partial 
wave and the phase difference between the 1 ^ and 2 ^ waves. The solid curves are the results of 
mass-dependent fits to the 7ri(1600) and 7r2(1670). 

In follow-up studies, COMPASS has looked at three-pion final states produced from a liquid 
hydrogen target. They have presented results on both 7 T~p P7 t~7t^7t~ and 7 T~p 158 . 

These studies confirmed the 1 ^ exotic signal found in the Pb running, but the production rates 
on Pb were found to be much higher than on hydrogen. Amplitude analysis showed that the exotic 
wave was produced with an = 1+ exchange in both reactions, but that the M — 1 in hydrogen 
is suppressed relative to the M = 0. Fig. ^ show the intensity of the 1 ^ wave in both the 

exotic signal in the 


TT TT ' TT 


and TT TT^TT^ systems. 


Fig. show the intensity of the 1 
Essentially the same information on the 1 ^ 
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Figure 27: (Color on line.) COMPASS results showing the intensity of the exotic 1 ^ wave. 
The solid (red) curve shows a fit to the corresponding resonances. The dashed (blue) curve is the 
7ri(1600) while the dotted (magenta) curve is background. (Figure reproduced with permission 
from reference (ill.) 



Figure 28: (Color on line.) COMPASS results showing the phase difference between the exotic 
1 ^ wave and the 2 ^ wave. The solid (red) curve shows a fit to the corresponding resonances. 
(Figure reproduced with permission from reference 153| .) 


(Stt) system can be found in references 160,162,163 


In order to more fully understand the three-pion system, COMPASS has studied their very 
large data set on three pions. In particular, the non-exotic 1++ partial wave which has historically 
been modeled using both the broad ai(1260) as well as non-resonant Deck-effect terms. They have 
reported on a new , nar ro w st ru ctur e in th is partial wave, the ai(1420), which decays dominantly 

. They report a mass of MeV and a width of 

this new state, results of the exotic 7ri(1600) have been 


to /o(980)7r 


153^^3 MeV 


164 


170 


165 


166 


167 


168 


Since their discovery o 


limited until a more detailed analysis can be completed. They report that an exotic signal in the 
1+ partial wave at 1.6 GeV/c^ is observed and that it shows a clean phase motion with respect to 
well-known resonances. Their results are consistent for both the and 7 t^7t~7t~ final states. 
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Figure 29: (Color on line) (left) The intensity = 1 ^ partial wave in the two ( 37 r“) final 
states. In the 1.6 GeV mass region, both the 7 r“ 7 r+ 7 r“ and the final states show the same 

structure. The peaks in the lower-mass region are unstable and probably artifacts of the analysis. 
The phase difference (right) between the and 1 ^ partial waves. (Figure reproduced with 
permission from reference [158| .) 


The state is observed in the pir decay mode for both charged and neutral ps. However they exclude 
that a narrow (150 MeV to 200 MeV width) state exists |166| , [169 , |170j . 

The COMPASS Collaboration has also examined the 7 T~p r]7T~p and 7 T~p p'7v~p reac¬ 
tions [134 , 159 . Amplitude analyses of both systems shows that the even partial waves are pro¬ 


duced with similar strengths in both system, while the odd partial waves are suppressed in the pn 
system relative to the pV system. The a2(1320) and a4(2040) are produced in both systems with 
a relative ratio given by phase-space factor and the pseudoscalar mixing angle. The odd partial 
waves, = (1,3, 5) ^ all correspond to exotic quantum numbers, and are produced much more 

strongly in the p' system. This is shown for the spin 1 and 3 partial waves in Fig. [^ In terms 
of quark structure, the p is predominantly an octet state, while the p' is predominantly a singlet 
state. Both VES and E852 reported that the 1 ^ exotic wave is the dominant feature in the p'lv 
final state, while it is strongly suppressed in the p channel. In addition, the relative phase motion 
between the exotic wave and the a 2 wave are different between the p and p' channels. Even though 
COMPASS saw the exotic 1 ^ wave in p'lr as the dominant wave, they were unable to confirm the 
resonant nature of the signal. Large non-resonant contributions appear to be required to describe 
both the intensity and the phase simultaneously. Essentially the same information on the 1 ^ 
exotic signal in the 77 system can be found in references 160 162 . 


5.3.7 Summary of 7ri(1600) results 


Table 21 summarizes the masses and widths found for the 7ri(1600) in the four decay modes and 
from the experiments which have seen a positive result. While the p'lr^ fiir and biir decay modes 
appear to be robust in the observation of a resonant 7ri(1600), there are concerns about the Sir 
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Figure 30: (Color on line) The intensity of the 77 ^“ (black) and r]7v~ (red) signals in the 

partial wave (a) and the 3 ^ partial wave (b). The signal for the t/tt” is much smaller than that 

for the 77 ' 7 r“ for all odd partial waves. (Figure reproduced with permission from reference [i^ .) 


final states. While we report these in the table, the results should be taken with some caution. 


Mode Mass (GeV) Width (GeV) Experiment Reference 


bi7T 

1.58 ±0.03 

0.30 ± 0.03 

VES 

171 

biTT 

1.61 ± 0.02 

0.290 ±0.03 

VES 

141 

biTT 

~ 1.6 

~ 0.33 

VES 

125 

biTV 

1.56 ±0.06 

0.34 ± 0.06 

VES 

126 

flTT 

1.64 ±0.03 

0.24 ±0.06 

VES 

126 

tj'tv 

1.58 ±0.03 

0.30 ± 0.03 

VES 

171 

r]'7T 

1.61 ± 0.02 

0.290 ±0.03 

VES 

141 

rf'iv 

1.56 ±0.06 

0.34 ± 0.06 

VES 

126 

piT 

1.593 ±0.08 

0.168 ± 0.020 

E852 

143 

tj'tv 

1.597 ±0.010 

0.340 ± 0.040 

E852 

145 

flTT 

1.709 ±0.024 

0.403 ± 0.080 

E852 

146 

bi7V 

1.664 ±0.008 

0.185 ±0.025 

E852 

147 

biTT 

~ 1.6 

~ 0.23 

GEAR 

149 

piT 

1.660 ± 0.010 

0.269 ±0.021 

COMPASS 

153 

tj'tv 

1.670 ±0.030 

0.240 ±0.050 

CLEO-c 

150 

all 

l.UUZ_g Qgg 

0.241 ± 0.040 

PDG 

118 


Table 21 : Reported masses and widths of the 7ri(1600) along with the 2014 PDG average. 


Models for hybrid decays predict rates for the decay of the tti. Equation]^ gives the predictions 
from reference Eg. 


Tihi : 7r/i : TTp : r]7T : itt]' 


170 : 60 : 5 - 20 : 0 - 10 : 0 - 10 

A second model from reference El predicted the following rates for a 7ri(1600). 


( 48 ) 
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Tvbi OTT 7r/i 77(1295)7r K*K 

PSS 24 9 5 2 0.8 

IKP 59 8 14 1 0.4 

These can be compared to the results from VES in equation [47| which are in moderate agreement. 
The real identification of the 7ri(1600) as a hybrid will almost certainly involve the identification 
of other members of the nonet: the r]i and/or the r][^ both of which are expected to have widths 
that are similar to the tti. For the case of the 771 , the most promising decay mode may be the fir] 
as it involves reasonably narrow daughters. 

We believe that the current data support the existence of a resonant 7ri(1600) which decays 
into 5 i 7 r, /itt and ryV, however, confirmation of the biir and /itt modes by COMPASS would be 
useful. For the pir decay, things are somewhat uncertain. As noted earlier, the phase motion results 
observed by both E852 and E852-IU are can be interpreted as either the 7r2(1670) absorbing the 
7ri(1600), or leakage from the 7r2(1670) generating a spurious signal in the 1 ^ channel. While 
COMPASS does not yet have final results, they seem to confirm that the exotic partial wave near 
1600 MeV does couple to the pir decay mode and they exclude a narrow-resonance interpretation 
of the 7ri(1600). The careful follow-on studies from COMPASS to more broadly explore the model 
space and production mechanisms have started to yield very interesting new results and we look 
forward to new information in the near future. 


5.4 The 7ri(2015) 

5.4.1 E852 Results on the 7ri(2015) 


The E852 experiment has also reported a third tti state seen decaying to both /itt fl46l and to 


biTT |147| . In the /itt final state, the 7ri(2015) is produced with = 1+ in conjunction with 
the 7ri(1600). The description of the 1 ^ partial wave requires two poles. They report a mass 
of 2.001 ± 0.030 ± 0.092 GeV and a width of 0.333 ± 0.052 ± 0.049 GeV. Fig. [31] shows the E852 
data from this final state. Parts e and / of this show the need for the two-pole solution. VES 
also examined the /itt final state, and their intensity of the 1 ^ partial wave above 1.9 GeV (see 
Fig. 18) is not inconsistent with that of E852 [126| . However, VES made no comment on this, nor 
have they claimed the existence of the 7ri(2015). 

In the bill final state, the 7ri(2015) is produced dominantly through natural parity exchange 
(M^ = 1+) while the 7ri(1600) was reported in both natural and unnatural parity exchange, where 
the unnatural exchange dominated. They observe a mass of 2.014 ±0.020 ±0.016 GeV and a width 
of 0.230 ±0.032 ±0.073 GeV, which are consistent with that observed in the /itt final state. Fig. [3^ 
shows the intensity distributions for several partial waves in this final states. The need for two 
states is most clearly seen in panel ( 6 ). VES also looked at the bin final state, but did not observe 
1 ^ intensity above 1.9 GeV fl 2 ^ . However, the intensity shown in Fig[T7|may be consistent with 
that observed by E852. The reported masses and widths are summarized in Table We note 
that this state does not appear in the summary tables of the PDG |118| . 


5.4.2 Interpretation of the 7ri(2015) 

With so little experimental evidence for this high-mass state, it is difficult to say much. We note 
that the observed decays, fin and bin are those expected for a hybrid meson. We also note 
that the production of this state is consistent (natural parity exchange) for both of the observed 
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Figure 31: (Color on line.) The fin invariant mass from E852 146| . (a) The partial wave 
(ai(1270)), (b) the 2 ^ partial wave (7r2(1670)) and (c) the exotic 1 ^ partial wave. The dotted 
(red) curves show the fits of Breit-Wigner distributions to the partial waves, (d) shows the phase 
difference between the 2 ^ and 1 ^ partial waves, while (e) shows the difference between the 

and 1 ^ partial waves. The dotted (red) curves show the results for a single tti state, the 7ri(1600). 
(f) shows the same phase difference as in (d), but the dotted (red) curve shows a fit with two poles 
in the 1 ^ partial wave, the 7ri(1600) and the 7ri(2015). (Figure reproduced with permission from 
reference fl46] .) 



Figure 32: The bin invariant mass from the E852 experiment, (a) shows the 1 ^ bin partial wave 
produced in natural parity exchange (M^ = 1^) while (b) shows the 1 ^ bin partial wave produced 

in unnatural parity exchange (M^ = 0“). In (c) is shown the 2++ ojp partial wave, while (d) shows 
the (jjp partial wave. The curves are fits to the 7ri(1600) and 7ri(2015) (a and b), the a2(1700) 
in (c) and the a4(2040) in (d). (Figure reproduced with permission from reference (MZl-) 


final states. In the case that the 7ri(1600) is associated with the lowest-mass hybrid state, one 
possible interpretation of the 7ri(2015) would be an excited state (as suggested by recent LQCD 
calculations [^). The mass splitting is typical of radial excitations observed in the normal mesons. 
In the case that the 7ri(1600) is identified as something else, the 7ri(2015) would be a prime candidate 
for the lightest mass hybrid. 
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Mode 

Mass (GeV) 

Width (GeV) 

Experiment 

Re 

hrence 

/iTT 

2.001 ± 0.030 

0.333 ± 0.052 

E852 


146 


biTT 

2.014 ±0.020 

0.230 ± 0.032 

E852 


147 



Table 22: Reported masses and widths of the 7ri(2015) as observed in the E852 experiment. The 
PDG does not report an average for this state. 

6 Outlook 

Hybrid mesons should be observed in nonets, and lattice QCD has recently made quite clear predic¬ 
tions that several of these nonets should have exotic quantum numbers. The lightest supermultiplet 
of hybrids, 0 1 1 ,2 contains one exotic quantum number nonet (1 ^). Excitations of 

these should contain another 1 ^ nonet, in addition to a single spin-0 (0^ ) and two spin-2 (2^ ) 

nonets. Experimentally, we have seen evidence for up to three isospin-1, spin-1 exotic states: 
7ri(1400), 7ri(1600) and 7ri(2105). As we noted it seems unlikely that the lightest of these is a 
hybrid, while the heavier two states could map onto the isospin-1 members of the two 1 ^ nonets 

predicted by lattice QCD. Unfortunately, lacking experimental evidence for other members of these 
nonets, it is difficult to draw solid conclusions. In particular, observation of one or both isospin-0 
members of either nonet is needed. Beyond the spin-1 exotic hybrids, observation of states with 
quantum numbers 2^ and or 0^ is needed. 

As we move forward over the next several years, new experimental efforts will join COMPASS in 
searching for these states. At Jefferson Lab, the ClueX experiment expects to take first physics data 
in late 2015. ClueX will use linearly polarized photons incident on a hydrogen target to produce 
these states. In the next five years or so, PANDA at FAIR will join these efforts with antiproton 
beams - hopefully dramatically new information will be gleaned from these experimental efforts. 

Keeping pace with experiment will require a number improvements in theory. It is hoped, for 
example, that lattice gauge theory will be able to provide predictions for light and charmonium 
hybrid masses on unquenched lattices with physical pion masses and multihadron operators rela¬ 
tively soon. Conclusive and comprehensive computations of strong and electromagnetic transitions 
of light and charmonium hybrids would also be most welcome. 

A new breed of models that is capable of reproducing central lattice results is also required. 
Ideally these models will reproduce the gluonic adiabatic potentials and the spectrum of heavy and 
light hybrids reasonably well. Presumably this will require a formalism that captures short range 
and long range dynamics in an approximate fashion without double counting or other conceptual 
issues. Perhaps a promising approach would be a model based on Coulomb gauge QCD with a vari¬ 
able number of constituent gluons and many-body interactions. Such a model should also be able 
to describe strong and electromagnetic decays reasonably accurately. These are demanding criteria, 
but new experimental and lattice data should provide many valuable clues, which will hopefully 
lead to a quantitative and qualitative understanding of this enigmatic sector of the Standard Model. 
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